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INTRODUCTION 

DPSI's  contract  with  the  Convective  Cloud  Physics  Branch 
(LYC) , Meteorology  Division  of  the  Air  Force  Geophysics  Labora- 
tory (AFGL)  F19628-76-C-0051  required  the  company  to  continue 
the  development  of  data  processing  techniques  which  are  em- 
ployed in  the  analysis  of  cloud  physics  data. 

This  report  is  a description  of  the  effort  expended 
by  DPSI  for  AFGL (LYC)  in  that  fulfillment.  We  have  organized 
our  work  in  sections  2 through  8 as  follows: 

Section  2:  Additions  and  alterations  using  the  PMS 
ID  System  as  a means  of  input 

Section  3:  Development  of  new  techniques  which 
employ  the  PMS  2D  System  as  input 

Section  4:  Development  of  new  techniques  which 
utilize  the  LEARJET  data  collection  system 

Section  5:  Alterations  in  the  various  radar  pro- 
grams which  have  been  necessitated  by  radar  hardware 
changes . 

Section  6:  Modifications  of  the  programs  for  the 
PDP8  on  board  computer 

Section  7:  Generation  of  new  analysis  programs 
for  the  benefit  of  LYC  scientists 


Section  8:  The  operating  instructions  for  the 
above  programs. 
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PMS  ID  PARTICLE  SIZING  SYSTEM 

The  one-dimensional  particle  sizing  system  on  board  each 
aircraft  collects  particle  size  distribution  data  for  one  sec- 
ond intervals.  Data  is  recorded  in  real  time  using  a Kennedy 
incremental  7 track  tape  recorder.  These  tapes  are  shipped 
to  AFGL/LYC  for  processing  and  analysis.  The  actual  details 
of  the  hardware  utilized  may  be  found  in  the  following  two 
reports;  "Development  and  Application  of  Data  Processing 
Techniques  and  Analytic  Procedures  to  Cloud  Physics  Data"  by 
DPSI,  and  "A  Manual  for  the  Particle  Sizing  Spectrometer  Sys- 
tem" by  PMS  Inc. 
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2.1  ID  PROCESSOR  (KNOLL1D) 

Program  KNOLL ID  processes  the  one-dimensional  particle 
data  and  calculates  parameters  such  as  particle  number  den- 
sity, liquid  water  content,  radar  reflectivity,  median  volume 
diameter,  and  many  others.  This  program  may  be  run  repeatedly 
with  different  particle  types  (ice  crystals)  and  different 
processing  options.  KNOLL ID  has  a full  edit  capability  which 
allows  the  counts  of  selected  channels  to  be  corrected  during 
specified  time  periods.  The  program  also  produces  various 
output  files  (mag  tape  and/or  cards)  to  be  used  with  additional 
processing  programs. 

This  program  is  actually  the  "workhorse"  of  LYC.  Since 
the  nature  of  all  the  work  performed  is  primarily  research, 

DPSI  is  often  required  to  make  alterations  in  order  to  ex- 
amine changes  in  the  scientist's  theories,  of  course,  modi- 

\ 

fications  are  also  required  from  time  to  time  because  of 
hardware  malfunctions  on  the  aircraft. 

During  the  past  year  ten  separate  versions  of  KNOLLlD 
have  been  maintained.  These  versions  were  required  to  meet 
the  specific  requirements  of  LYC  personnel.  Details  of  the 
basic  program  (v2.00)  are  described  in  a separate  report  en- 
titled "KNOLLlD  OPERATORS  MANUAL".  This  report  explains  the 
j different  operating  versions,  and  the  required  control  card 

changes  for  each. 

Version  2.01  was  written  to  correct  for  the  upper  alti- 
tude pressures  of  flight  A75-19  (15  APR  75);  it  should  be 
used  for  this  flight  only.  The  modifications  include  the 
following : 
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!•  Pcor  " 1-25‘P-llO 


where 


P = corrected  pressure 
cor 

= calculated  pressure 


2.  For  aircraft  time  between  12:07:08  and 
12:48:53  use 

Ap  = 37.06  mb  FIXED 

Version  2.02  was  written  to  allow  processing  of  data 
with  a faulty  pressure  sensor.  This  version  will  accept  a 
time-height  profile  as  part  of  the  input  deck.  The  height 
at  a given  time  is  determined  by  linear  interpolation. 

The  pressure  is  then  calculated  from  the  interpolated  height 
using  the  following  equation: 


P = P0[l  - (^-)H T 1/aR| 


p0'  T0/  R and  a are  also  input  variables  and  are  defined  as 
follows : 

Pq  = Pressure  at  ground 
T0  = Temperature  at  ground 
R = Gas  constant 
a = Lapse  rate  at  ground 
HT  = Height 

Version  2.03  was  written  at  the  request  of  Dr.  Robert 
M.  Cunningham  of  the  Meteorology  Division,  AFGL.  These 
changes  were  experimental  in  nature  and  intended  for  flight 
A75-19  only.  This  version  changed  the  type  2 (wet  snow) 
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breakpoint  to  5 instead  of  3 for  the  precip  probe  data. 

Refer  to  the  KN0LL1D  OPERATORS  MANUAL  for  a complete  des- 
cription of  the  breakpoint  method  used  for  the  equivalent 
melted  diameter  determination.  Similar  cloud  probe  changes 
were  accomplished  with  data  cards. 

Versions  2.04  and  2.05  were  written  simultaneously  and 
include  a complete  data  modification  routine  as  requested  by 
Dr.  Robert  M.  Cunningham.  V2.04  was  necessitated  to  apply 
these  modifications  to  the  A75-19  flight.  V2.05  is  used  for 
all  standard  flights.  The  goal  of  these  modifications  was, 
ultimately  to  increase  the  liquid  water  content  and  radar 
reflectivity  because  of  certain  hardware  design  shortcomings. 
The  details  of  the  techniques  utilized  along  with  examples 
will  be  found  in  the  KNOLL ID  OPERATORS  MANUAL. 

The  justification  for  this  modification  is  twofold. 
Firstly,  there  are  "blind  spots"  associated  with  the  ID  in- 
strument. A blind  spot  may  be  defined  as  that  portion  of 
the  spectra  which  lies  between  the  limits  of  two  probes. 

When  the  instrument  is  measuring  water  droplets  there 
is  no  blind  spot  present.  However,  it  is  the  ice  crystal 
shadow  length,  being  reduced  to  equivalent  melted  diameter, 
that  causes  the  gap.  These  blind  spots  are  a function  of 
particle  type. 

For  example,  when  measuring  needles  the  upper  limit  of 
the  cloud  probe  is  142  microns  and  the  lower  limit  of  the 
precip  probe  is  241  microns,  then  the  ID  instrument  is  "blind" 
to  particles  in  the  142  to  241  micron  range.  Obviously  this 
could  have  a serious  effect  on  water  content  and  reflectivity 
calcul utions . 


9 


Secondly,  there  is  a design  feature  of  the  ID  system 
which  requires  modification  of  the  data  sampled.  This  fea- 
ture is  known  as  "end-rejection".  The  reason  for  this  re- 
jection is  the  philosophy  of  the  ID  system;  if  the  ending 
diode  is  occluded,  there  is  no  way  to  estimate  the  true 
particle  length  and  it  was  felt,  at  the  engineering  design 
level,  that  it  would  be  better  to  eliminate  the  particle 
rather  than  counting  it  as  one  with  a lesser  diameter.  With 
that  in  mind,  the  end  diode  rejection  feature  was  incorporated 
into  the  ID  system.  An  area  of  concern,  however,  is  the  num- 
ber of  particles  being  rejected;  with  the  present  ID  system 
there  is  no  way  to  determine  this.  It  may  be  of  consider- 
able consequence  because  the  larger  particles  have  the  great- 
est probability  of  being  rejected,  and  it  is  precisely  these  • 
particles  which  will  contribute  heavily  to  the  liquid  water 
content  and  radar  reflectivity. 

Since  the  particle  spectra  is  assumed  to  be  continuous, 
a particle  count  of  zero  in  a given  channel,  flanked  by  non 
zero  channel  counts  would  indicate  the  zero  channel  value  is 
incorrect.  There  is  a high  probability  that  there  really 
were  particles  in  the  zero  channel,  but  they  were  rejected 
because  the  ending  diode  was  occluded.  This  result  will  also 
seriously  affect  the  water  content  and  reflectivity  calcula- 
tions. 

These  versions  (2.04  and  2.05)  also  allow  for  nine  dif- 
ferent particle  types  to  be  processed.  Two  of  these  are  ex- 
perimental and  are  used  exclusively  by  Dr.  Robert  M.  Cunning- 
ham. These  types  are  designated  as  RMC1  and  RMC2.  The  com- 
plete list  includes  the  following: 


RAIN 


WET  SNOW 


LARGE  SNOW 


SMALL  SNOW 


BULLET-ROSETTES 


COLUMNS  (4:1) 
NEEDLES  (7.5:1) 
RMC1 


RMC2 


Note  only  odd  particle  type  numbers  are  specified.  All  even 
numbers  are  used  internally  within  the  program. 


The  equivalent  melted  diameter  calculation  has  also  been 
modified.  This  method  allows  two  breakpoints  to  be  specified 
for  each  particle  type.  The  first  breakpoint  is  specified  in 
units  of  N (channel  number) . This  controls  the  N'  equations: 


adjusted  channel  number 

channel  number 

particle  type  code 

channel  number  breakpoint 

intercept 

slope 


The  crystal  size  is  then  calculated  using 


crystal  size  (mm) 
probe  diode  width  (mm) 
adjusted  channel  number 


After  the  crystal  size  has  been  determined  the  ecuivalen 
melted  diameter  is  calculated  for  each  channel  using  the 
equations : 


where 


equivalent  melted  diameter  (mm) 
crystal  size  (mm) 
crystal  size  breakpoint 
particle  tvpe  code 


exponen 


bote  that  steps  1 and  3 actually  allow  for  two  ecuations. 
The  equation  chosen  is  de"endent  upon  the  channel  number 
(in  step  1)  or  the  crystal  size  (in  step  3).  Refer  to  the 
KI-iOLL ID  OPERATORS  MANUAL  for  all  the  coefficients  and  ex- 
ponents used  in  these  calculations . Also  included  are  the 
channel  number  (class),  adjusted  channel  number,  crystal 
size  (SC  for  scatter,  CL  for  cloud  and  PP.  for  rroepr)  , and 
equivalent  melted  di  ameter  for  each  size  clas^  . The  second 
breakpoint  is  specified  in  units  of  mm  (shadow  length). 


This  controls  the  equivalent  melted  diameter  equation  as 
shown  in  step  3.  Refer  to  the  KN0LL1D  OPERATORS  MANUAL  for 
a deta iled  discussion  of  both  particle  typing  and  equivalent 
melted  diameter  calculation. 

An  additional  parameter  is  also  calculated  in  these 
versions.  Kappa  (or  K)  is  defined  as  M//z.  K is  calculated 
four  times  per  averaging  interval;  once  for  each  probe  plus 
the  cloud  and  precip  probes  combined.  In  each  case  the  aver- 
age mass  and  average  radar  reflectivity  are  used.  The  plot 
tape  produced  (input  to  KNPLT1D)  has  been  changed  to  include 
these  four  values. 

In  addition,  to  conserve  space  and  reduce  calculation 
time,  the  six  moments  of  distribution  previously  calculated 
have  been  eliminated.  The  algorithm  used  in  these  calcula- 
tions has  been  retained  for  future  use,  if  necessary. 

The  RAPP  tape  produced  (input  to  RAPP)  has  also  been 
changed  in  these  versions.  The  new  format  is  considerably 
shorter  and  allows  for  more  efficient  processing.  The  new 
format  contains  9 words  pet  record.  Items  included  on  the 
new  tape  are:  date,  time,  averaging  interval,  liquid  water 
content  (for  each  probe) , and  radar  reflectivity  (for  each 
probe) . 

During  December  1975  certain  hardware  modifications 
necessitated  a change  in  the  Kennedy  tape  format.  This 
new  format  is  summarized  on  the  following  page. 


Word  # 

NEW 

OLD 

1 

elapsed  time 

sync 

17 

status 

word 

( scatter ) 

probe  configuration 

33 

status 

word 

(cloud) 

elapsed  time 

49 

status 

word 

(precip) 

unused 

A hardware  change  in  the  Knollenberg  ID  system  required 
changing  the  location  of  the  elapsed  time  word  within  the 
tape  record.  In  addition,  the  high  order  digit  has  been  lost 
A routine  was  added  to  this  version,  to  recreate  that  high 
order  digit. 

Vei'sion  2.06  was  written  to  accommodate  these  changes. 
This  became  the  standard  operating  version  for  all  flights 
after  December  1975. 

At  the  request  of  Dr.  Robert  M.  Cunningham  the  data 
modification  routine,  used  until  this  time,  was  modified 
and  incorporated  into  versions  2.07  and  2.08.  The  changes 
were  primarily  initiated  when  the  data  of  flight  A76-04 
(11  FEB  76)  was  processed.  Essentially  the  change  elimi- 
nated the  entire  data  modification  routine  if  there  is  in- 
sufficient data  from  the  precip  probe  (i.e.  less  than  two 
channels).  Version  2.07  became  the  standard  (from  2.06) 
and  Version  2.08  was  created  from  2.02.  This  allowed  the 
EDP  runs  (with  the  time-height  profile)  to  be  processed  with 
the  revised  data  modification  routine. 

Version  2.09  includes  a line  printer  nlot  of  median 
volume  diameter,  liquid  water  content,  and  radar  reflec- 
tivity as  a function  of  time  for  each  pass.  The  parameters 
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are  total  values,  i.e.  cloud  and  precip  probe  combined; 
plotted  on  a logarithmic  axis  81  columns  long.  A range  of 
+ 4 is  acceptable  with  a resolution  of  0.1.  Data  outside 
these  limits  is  plotted  with  a special  character  either 
above  or  below  the  plot. 

Tiie  flowchart  on  the  following  page  depicts  the  "evo- 
lution" of  the  ten  versions  maintained  during  the  past  year. 
These  programs  are  permanently  archived  (with  backup)  and 
may  be  retrieved  for  processing  or  modification  at  any  time. 
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2.2  ID  PLOTTER  (KNPLT1D) 


The  plotting  program  for  the  PHS-lD  device  has  undergone 
numerous  changes  during  the  past  year.  Six  different  types 
of  plots  are  produced  displaying  the  processed  data  in  dif- 
ferent forms.  The  input  tape  used,  produced  by  program 
KN0LL1D,  contains  averaged  data,  one  record  per  average  in- 
terval. In  addition  the  program  can  read  a processed  tape 
containing  Learjet  data. 

Each  plot  tape  to  be  produced  requires  one  plot  request 
card.  This  card  contains  the  information  necessary  to  pro- 
duce the  plot  (start  time,  stop  time  etc.)  with  the  appro- 
priate title.  The  details  of  the  plot  request  cards  are 
included  in  the  operating  instructions. 

To  reduce  plotting  time  and  paper  consumption  the  plots 
are  produced  on  35  mm  film.  The  following  sections  explain 
in  detail  the  plots  available  and  the  calculations  performed. 
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2.2.1  Z-li  scatter  diagrams 

This  option  produces  two  scatter  diagrams  of  mass  (M) 
in  gms/X3  vs.  reflectivity  (Z)  in  nm6/m3  on  a full  log 
scale.  The  first  plot  used  precip  probe  data  only,  the 
second  uses  cloud  and  precip  combined.  Any  mass  value  out- 
side the  range  of  10~4  through  10+1  is  not  plotted.  The 
reflectivity  limits  are  10-4  and  10+5. 

A least  square  logarithmic  fit  is  calculated  and  drawn 
through  the  data.  The  equation  of  this  line  is  included  on 
the  plot  in  the  form 

M = a Zb 

The  output  listing  for  this  plot  type  includes  the 
equation  coefficients  (both  linear  and  exponential),  the 
average  Z and  14  values,  with  their  standard  deviations. 
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2.2.2  Z— 1-1  histograms 

This  option  produces  two  histograms;  the  first  mass  (M) 
vs  time  and  the  second  reflectivity  (Z)  vs  time.  The  plots 
are  produced  using  a semi-log  scale.  The  Z and  M plot 
limits  are  the  sane  as  those  defined  in  section  .2.1.  The 
time  axis  is  expressed  in  units  of  seconds  from  a given 
start  time.  These  plots  restrict  the  data  to  five  minutes 
only.  Any  data  exceeding  this  maximum  will  be  ignored. 

Each  plot  has  three  traces,  one  per  probe,  with  a dif- 
ferent plotting  symbol  used  for  each  probe.  A square  indi- 
cates data  from  the  scatter  probe,  a circle  represents 
cloud  data,  and  a triangle  is  used  for  precip  data.  The 
output  listing  shows  only  the  number  of  points  plotted  per 
trace . 


J 

i 
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2.2.3  VCO  histograms 

This  option  produces  three  VCO  histograms,  each  con- 
sisting of  two  lines.  The  VCO's  plotted  are:  pressure, 
temperature,  heading,  dev/point,  acceleration,  and  J-W  water 
content.  The  plots  are  used  exclusively  for  system  veri- 
fication and  are  not  published  outside  the  laboratory.  For 
this  reason  there  is  no  title,  pass  or  other  identification 
on  the  plot.  The  time  axis  is  similar  to  the  Z-M  histograms 
and  also  has  the  same  five  minute  maximum. 

The  acceleration  and  water  content  plots  have  fixed 
scales;  the  others  have  a sliding  scale  with  a fixed  maxi- 
mum range.  The  limits  used  for  the  sliding  scale  plots 
are  a function  of  the  data  to  be  plotted.  They  are  deter- 
mined automatically  by  the  program.  Figure  2 below  shows 
the  pertinent  axis  information.  The  output  listing  for 
this  plot  shows  only  the  number  of  points  plotted. 


PLOT# 

VCO 

UNITS 

SYMBOL  ' 

AXIS  TYPE 

! AXIS  LIMITS 

1 

pressure 

acceleration 

mb 

g 

A 

- 

sliding 

fixed 

30  mb  max  | 

-1  to  +1 

h 

9 

temperature 

deg  C 

A 

sliding 

12°  max 

dewpoint 

deg  C 

+ j 

sliding 

12°  max 

heading 

deg 

A 

si iding 

60°  max 

! j 

JW-  LVJC 

gm/M  3 

fixed 

- . 1 to  + . 8 

Figure  2:  VCO  Histogram  Axis  Information 
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2.2.4  Density  spectra 


'1 

I 

' I 

r. 


•! 

;? 


This  option  produces  for  a selected  probe  six  plots 
per  page  of  the  log10  of  number  density  vs  the  equivalent 
melted  diameter.  The  six  pictures,  rather  than  one,  results 
from  the  fact  that  it  is  desired  to  examine  this  plot  sepa- 
rately as  liquid  water  content  ( LWC ) increases.  Thus  the 
procedure  is  to  first  determine  the  combined  LWC  for  the 
cloud  and  precip  probes  and  determine,  as  a function  of  the 
LWC,  which  of  the  six  plots  is  to  be  used  to  display  the 
data;  then  the  appropriate  log  15  of  the  number  density  and 
equivalent  melted  diameter  is  retrieved  to  produce  the 
plotting  points  within  the  plot  selected.  The  LWC  limits 
for  the  six  plots  are  defined  on  the  plot  request  card. 


After  all  this  data  has  been  plotted  additional  calcu- 
lations are  performed.  These  calculations,  (different  types 
of  averaging)  requested  by  Dr.  Robert  M.  Cunningham,  are 
shown  in  the  following  ten  steps. 

1.  Two  parameters  must  be  calculated  before  the  averages  can 
be  calculated. 

These  are  . . . 

ICNTX(N)  and 

SUMX(N)  for  N = classes  1 to  15 


ICNTX  (II) 


where 


SUMX(N)=  l 
i=l 


DENSITY  (N)  j_ 


ICNTX  (N)  = the  number  of  sarnies  in  each  class  N 


: ■ 


note:  The  normalized  density  is  being  used  at  this  time 

also  1CNTKX  is  the  maximum  ICNTX(N) 

2.  "INSTRUMENT  LOG  (AVE  DENSITY)" 

Calculated  for  each  non-zero  class  i.e.  only  if  ICNTX(N)>0 

AVE I NS  ( N ) = SUMX ( N ) / IC NTMX 

AVEINSL(N)  = LOG  (A VEINS  (N)  ) 

3.  "EXPONENTIAL  LOG  (AVE  DENSITY)" 

a.  Before  this  is  calculated  the  True  Ave  Log  Density  is 
determined  for  each  class  where  . . . 

ICNTX (N)  > . 9 ( ICNTMX) 

note:  JJPTS  = number  of  classes  where  above  relation  is  true 

using  the  equations 

A VETRU (N ) = SUMX (N) / ICNTX  ( N ) 

AVETRCJL  (N)  = LOG  (AVETRU  (N)  ) 

b.  The  set  of  data  points  (DIAM (N) , AVETRUL (N) ) 
for  ICNTX (N)  > 0 are  linearly  curvefit 

note:  IIPTS  = number  of  classes  where  above  relation  is  true 

c.  The  slope  (m)  and  intercept  (b)  are  printed  and  ex- 
pressed exponentially  to  satisfy  the  equation  ... 


N/M3-mm  = N e' 
o 


where  N 


2.3026 (b) 


and  X = 2.3026 (m) 

d.  The  exponential  Log  (AVE  DENSITY)  j.s  calculated  and 
plotted  for  the  full  spectrum  (classes  1-15) 

using  . . . 


AVEEXP(N)  = N e 
o 


> *D IAM  (N) 


AVEEXPL (N)  = LOG (AVEEXP (N) ) 


where  DIAM(N)  = center  diameter  of  class  N 


AVEEXPL  (N)  is  printed  only  for  classes  v/here  ICNTX(N)  > 0. 
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4.  "INSTRUMENT  MASS  AND  REFLECTIVITY" 


I 


a . INSTRUMENT  MASS 

Calculated  individually  for  each  class  where  ICNTX(N)>0 
using  the  instrument  average  density  (AVETKS) 

MINST (N)  = £ BW(N) *AVEINS (N) *DIAM(N) 3 

i 

where  BW(N)  = barwidth  of  class  N 

Also  the  cumulative  instrument  mass  is  calculated  as 

IIPTS 

TMI  = [ MINST (N) 

N=1 

b.  INSTRUMENT  REFLECTIVITY 

Calculated  individually  for  each  class  v/here  ICNTX(N)>0 
using  the  instrument  average  density  (AVEINS) 

Z INST (N)  = BW(N) *AVEINS (N) *DIAM(N) 6 


Also  the  cumulative  instrument  reflectivity  is  calculated  as 

IIPTS 

TZI  = l ZINST(N) 

N=1 


9 

| 

i 


5.  "EXPONENTIAL  .‘IASS  AND  REFLECTIVITY" 

a.  EXPONENTIAL  MASS 

Calculated  individually  for  each  class  where  ICNTX(N)>0 
using  the  exponential  average  density  (AVF.EXP) 

riEXP  (N)  = g-BV;  (N)  *AVE EXP  (N)  *DIAH  (N)  3 

Also  the  cumulative  instrument  mass  is  calculated  as 

IIPTS 

TME  - l EXP(M) 

N=1 

b.  EXPONENTIAL  REFLECTIVITY 

Calculated  individually  for  each  class  where  ICNTX(M)>0 
using  the  exponential  average  density  (AVEEXP) 

Z EXP ( N ) - BW ( N ) * AVEEX P ( N ) *DIAK ( N ) 6 

Also  the  cumulative  instrument  reflectivity  is  calculated  as 

IIPTS 

TZE  = l ZEXP(N) 

N’-l 


. Mi 
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6.  "CUMULATIVE  PERCENT  "ASS" 


CUMASS (M)  = MINST  (U)/TMI+TMASS(N-1) 


wnere 


N-l 

(TMASS (N-l) = l CUMASS (i) 
i=l 


7.  "TOTAL  MASS" 

To  account  for  integration  from  LM  to  ™ rather  than  from 
0 to  ”,  the  equation  from  SAMS* #2  p 58  is  modified. 


MTT  = M 


N / D 3 e ADdD 


( >D j ) - S 1 o 


'1 


MTT  = N e XDl 
6 o 


A 


3 3D  i 2 6Dj  6 

+ + + - 

A 3 A3 


A1*  I 


letting 


NTT  - 


N e 
o 


-ADj 


, since  the  slope  A should  always 
be  positive 


the  final  equation  becomes 


MTT  = J NTT 
6 


+ 


D j 6 


If  D]  =0  it.  can  be  seen  that  the  equation  is  icer.tical  to 
that  from.  SATIS**  2 p 58 


o 


A ! 


*Hydrom  r Pa ram  . Th  i 1 * or 

thi  SAMS  Rail  Erf  : n Prog  ; Plank,  v.  ' . , i . 


8 . "TOTAL  REFLEC TI\ l . * 


To  account  for  integration  from  D:  to  * rather  than  from 
0 to  <*>,  the  equation  from  SAJIS*  #2  p 53  is  modified. 

X 

2TT  = Z.  _ . = N / D6e~ 'DdD 

: D, c 6D i 5 30D  j 4 120D]  3 

2TT  = N e"'U>  + — + + + 

X X2  X3  X4 


360D’  2 720D)  720~j 


letting  N e-ADj 

NTT  = — 


the  final  equation  becomes 

6D 1 5 3 OD i u 1 2 0D i 3 360D i 2 720D,  7 


IIIUUP^  -ny" •£— ' 


21 


1C.  When  the  particle  type  being  analyzed  is  rain  only  the 
total  mass  and  total  reflectivity  should  be  integrated 
from  Dj  to  5 mm  only.  This  is  accomplished  by  sub- 
tracting the  5 mm  to  infinity  value  from  the  Dj  to 
infinity  value  for  M and  Z. 

a.  "MASS  CONSIDERATION" 

MTT  = J-  Nq  / D3e_ADdD  = £ Nq  { J D3e~XDdD -/  D3e~XDdD} 
Dj  Dj  5 


from  step  7 the  Dj  to  <*  evaluation  is 


30^  6D  j 6 


+ — ] = 


> 3 X4 


'XDl  [loom] 


using  the  same  analysis  the  5 nm  to  °°  evaluation  is 


'JA  [±-  + 


llll  + ll?-  + 1_]  =;  C"5A  [UPDM] 

X X 2 X 3 X 4 


and  finally 


MTT(Dl:5)  = Z No  f e"‘Ul (L0DM)-e'5X (UPDM)  } 


b.  "REFLECTIVITY  CONSIDERATION " 

Repeating  the  previous  analysis  for  reflectivity  the 
final  equation  becomes 


- 


29 


2.2.5  Median  volume  diameter  plots 

Tnis  option  produces  three  different  sets  of  scatter 
type  plots.  A set  consists  of  two  plots,  each  done  the  same 
way  only  using  different  data;  the  first  plot  uses  precip 
data  only  and  the  second  uses  precip  data  combined  with 
cloud  data.  The  three  sets  are  described  in  the  following 
sections . 

2. 2. 5.1  DQ  vs  (Z/M)1/3 

For  each  averaging  interval  this  plot  calculates  the 
ratio  of  reflectivity  to  mass,  and  then  takes  its  cube  root. 
This  value  is  checked  to  insure  it  is  within  the  horizontal 
axis  limits  of  10-1  and  1Q+1 . The  median  volume  diameter, 
for  the  same  interval,  is  checked  to  insure  it  is  within  the 
vertical  axis  limits  of  10“ 1 and  10+1 . 

After  the  data  is  plotted,  a least  square  logarithmic 
fit  is  calculated  and  drawn  through  the  data.  The  equation 
of  this  line  is  included  on  the  plot  in  the  form 

Dq  = a((Z/M)1/3)b 

The  output  listing  for  this  plot  includes  the  equation 
coefficients  (both  linear  and  exponential) , the  average 
( Z/M) 3 and  Dq  values,  and  their  standard  deviations. 


This  plot  is  similar  to  the  one  previously  discussed 
the  only  difference  being  the  x-variable  calculation.  In 
this  case  the  ratio  of  mass  to  square  root  of  reflectivity 
is  calculated.  The  remaining  calculatons  axis  scales,  and 
output  are  the  same  as  the  previous  plot. 

2. 2. 5. 3 ND  VM  vs  D/D 
o'  o 


The  calculations  performed  after  the  plot  is  produced 
were  requested  by  Dr.  Robert  M.  Cunningham.  The  procedure 
is  a "follow-up"  to  the  data  modification  technique  shown 
in  section  7.1.10.  This  plot  differs  primarily  from  the 
previous  two  in  that  it  uses  every  non-zero  channel  within 
an  averaging  interval.  That  is,  for  each  averaging  interval 
as  many  as  15  points  may  be  plotted.  (If  cloud  and  precip 


data  are  being  used,  then  there  may  be  as  many  as  30  points 
plotted  per  averaging  interval.) 


The  points  are  plotted  in  the  following  manner.  For 
each  averaging  interval  the  15  (or  30)  class  diameters  are 
divided  by  the  median  volume  diameter.  These  become  the 
x-coordinates  of  the  plotted  points.  The  y-coordinate  is 
determined  by  multiplying  the  number  density  of  each  class 
times  the  fourth  power  of  the  median  volume  diameter.  This 
product  is  then  divided  by  the  calculated  mass  of  the  inter- 
val. If  the  density  is  zero,  there  is,  of  course,  no  point 
plotted.  The  plot  is  produced  on  a semi-log  grid,  with  the 
x-axis  being  linear.  After  all  the  points  are  plotted,  a 
least  square  exponential  fit  is  calculated  and  drawn  through 
tne  data.  Tne  coefficients  of  the  equation  are  shown  on 


I 

* 


the  plot  in  the  form 


ND  4/M  = a eb  D/Do 


After  the  plot  is  produced,  the  calculations  performed 
are  outlined  below. 


1.  Divide  the  x-axis  into  21  banc’s 


2.  Each  point  plotted  is  categorized  by  its  x-coordinate 
into  one  of  these  bands. 

a.  The  y-coordinate  is  accumulated  with  other  values  in 
the  same  band . 

b.  The  number  of  points  in  each  banded  is  counted. 

c.  When  all  the  data  has  been  categorized  and  summarized, 
the  average  y value  for  each  band  is  determined. 

3.  The  bands  with  a zero  average  y value  are  examined.  This 
zero  results  because  of  one  of  two  conditions. 

a.  It  may  be  that  a particular  band  had  no  points  (these 
are  called  "uncounted"  or  type  1 zeroes) . 

b.  It  may  be  that  all  the  y-coordinate  for  a particular 
band  were  zero  (these  are  called  "counted"  or  type  2 
zeroes) . 


4.  Each  type  of  zero  may  appear  in  two  places: 

a.  as  a surrounded  zero  (or  consecutive  zeroes)  where 
they  are  surrounded  by  non-zero  bands  on  both  sides. 

b.  as  an  ending  zero  (or  consecutive  zeroes) 
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The  following  examples  show  both  cases. 


example  A.  surrounded  zero 


band# 


average  y 


1.24 


3.94 


1.85 


surrounded  zero 


surrounded  zeroes 


example  B. 


ending  zero 


bands 


average 


ending  zero 


3.28 


1.29 


ending  zeroes 


5.  An  attempt  is  made  to  eliminate  all  the  zeroes.  The 
elimination  technique  used  is  dependent  upon  the  type  (1  or  2) 
and  place  (ending  or  surrounded)  of  zero.  The  hierarchy 
and  tecnnique  used  to  eliminate  these  zeroes  is  shown  on  the 
following  page. 


log  Y log 


I 
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zero (type  and  place)  technique 


1) 

#2 

surrounded 

3 point  running  log  mean 

2) 

#2 

ending 

extrapolation 

3) 

#1 

ending 

extrapolation 

4) 

#1 

surrounded 

interpolation 

a . Type  2 surrounded  zeroes  are  removed  using  a three 
point  running  log  mean;  this  changes  the  surrounding 
non-zero  points  also. 


4 - 

X 

X 
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Figure  3 : Examples  of  type  2 surrounded  zero  elimination 
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b.  Ending  zeroes  of  types  1 and  2 are  eliminated  by  extra- 
polating from  the  equation 

y = ce[ax2  + bx] 

The  equation  is  a least  square  fit  of  four  non-zero 
points  immediately  preceding  or  following  the  zeroes 
to  be  eliminated. 


o = zero  point  x = non-zero  point 


the  fitting  equation 
uses  points  y3  y^ 

r f 

y6  and  y8  to  elimi- 
nate the  zeroes  in 
bands  1 and  2 


Figure  4 : Example  of  ending  zero  elimination 


Zeroes  at  the  other  end  are  eliminated  exactly  the 
same  way  using  the  four  non-zero  points  preceding 
the  zeroes  to  be  eliminated. 


c.  Type  1 surrounded  zeroes  are  removed  using  linear 
interpolation  betv/een  adjacent  non-zero  values. 


i 

i. 

v 


I 


; 


2 3 

band  # 


-6 


2 3 4 

band  # 


y 4-y2 

y3'  “ 72  + (^7}  (X3“X2> 


y'l  = 2 + (3-2 


y$  = 2 - h = ik 


° - zero  point  x = non-zero  point 


Y5-Y2 

y3  " Y2  + (X3'X2) 


1-3, 


y3'  = 3 + (£|)  (3-2) 


y3  = 3 - (|)1  = 2 1/3 


yC 


y 2 + 


y5-y2 


X r-x 


-)  (x4-  x2 ) 


yC  = 3 - (j)2  = 1 2/3 


Figure  5 : Examples  of  type  1 surrounded  zero  elimination 


6.  After  the  zeroes  have  been  eliminated  the  results  are 
shown  in  the  output  listing.  The  bands  that  had  a zero 
average  eliminated  are  flagged  with  a 1 or  2 indicating 
the  type  zero  removed. 
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2.2.6  VCO  plots 

This  option  will  plot  any  VCO  versus  any  other  VCO. 

With  this  option  only,  the  calculated  liquid  water  content 
(M)  and  calculated  reflectivity  (Z)  for  each  probe  are  con- 
sidered VCOs.  This  allows  selected  Z vs  height  or  M vs 
height  profiles  to  be  plotted.  The  user  has  complete  free- 
dom to  specify  any  VCO  on  either  axis.  Plots  may  be  either 
in  scatter  mode  or  line  mode,  at  the  user's  option. 

If  the  plot  of  JW-LWC  vs  HEIGHT  is  specified,  an  ad- 
justment option  is  available.  This  adjustment  allows  for 
all  the  data  points  to  be  shifted  by  some  reference  equation. 
To  invoke  this  option  certain  adjustment  parameters  are 
necessary.  These  parameters  are: 

L = number  of  levels 
XJ  = origin  of  level 
SL  = slope  of  level 
HT  = height  (meters) 

This  information  is  input  via  a namelist  card  as  shown  in 
the  operatinq  instructions . The  adjustment  performed 
follows . 


Height 


XJj  XJ2  XJ3  XJU 

JW-LWC 


Figure  6 : JW-LWC  vs  height  adjustment 

Determine  the  level  (I)  of  each  data  point  (J,H)  such 
that  HTj+^  > F > HTj . Calculate  the  adjusted  water  content 

(JW)  using  the  equation 

JV7  = J-fXJj  + SL  (HTj-H)  ] 

Note:  if  height  (H)  is  greater  than  the  maximum  height  (HTj ) , 
than  HTj  is  used  in  the  ecuation  for  HTt.  The  sign  reversal 
will  account  for  the  correct  JW  adjustment. 
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2.3  PMS-lD  UTILITY  PROGRAM 


KNIUTIL  is  a generalized  tape  utility  program  designed 
to  operate  on  the  PMS-1D  input  tape  in  order  to  verify  the 
operation  of  the  Knollenberg  PMS  ID  device.  There  are  five 
different  options  the  user  may  select.  These  options,  are 
listed  here  and  further  explained  in  the  operating  instruc- 
tions . 


1.  produce  a decimal  dump  of  a given  number  of  files; 

2.  same  as  1,  above,  except  that  the  dump  is  in  octal; 

3.  produce  a selective  decimal  dump  specifying  a 
particular  set  or  sets  of  records; 

4.  copy  a tape  specifying  files  and  records  to  be 
copied; 

5.  produce  a selective  probe  dump  which  will  list  on]y 
the  channel  counts  of  a selected  probe  or  the  VCOs. 


In  addition  to  these  options,  three  status  files  are 
created  which  may  be  listed  at  the  users  option.  These  files 
list  in  an  easily  readable  format  the  three  status  words  that 
appear  in  each  record. 
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3.  PMS  2D  PARTICLE  DISPLAY  SYSTEM 

The  two-dimensional  (2D)  particle  display  system  on  board 
each  aircraft  has  some  important  advantages  over  the  earlier 
ID  models.  Firstly  there  are  32  sensors  in  each  probe  re- 
cording particle  image  data.  The  ID  system  has  24  diodes 
with  a maximum  particle  size  capability  of  15  classes.  The 
cloud  probe  covers  the  size  range  of  25  to  800  microns;  the 
precip  probe  covers  200  to  6400  microns.  The  second  dimen- 
sion is  achieved  by  taking  readings  over  time  so  that  a two 
dimensional  picture  of  the  shadow  is  made.  The  sampling 
rate  is  adjusted  to  the  speed  of  the  aircraft  so  that  a read- 
ing would  be  taken  every  25y  of  length.  That  is,  if  the  air- 
craft flies  at  100  meters/sec,  the  sampling  rate  would  have 
to  be  4 megahertz.  This  exact  ratio  can  not  be  maintained 
perfectly,  so  the  results  are  modified  slightly  in  the  com- 
puter according  to  the  true  air  speed  of  the  aircraft.  Like 
the  ID  device,  the  output  is  turned  on  when  a sensor  is  shut 
off,  and  continues  until  all  the  sensors  are  back  on,  but 
this  device  will  output  the  status  of  each  of  the  32  sensors 
every  four -millionth  of  a second*  until  all  the  sensors  are 
back  on.  Thus  the  2D  device  gives  a picture  of  the  particle (s) 
as  subsequent  readouts  are  place  together,  and  will  not  give 
incorrect  results  when  two  particles  are  seen  simultaneously. 

An  additional  advantage  of  the  2D  system  is  the  elimi- 
nation of  the  end-rejection  technique  utilized  in  the  ID 
system.  If  a particle  occludes  either  ending  diode,  it  is 
still  recorded.  These  particles  may  then  be  included  or 
rejected  manually  with  the  appropriate  software. 

* 


I 


i 


for  the  Learjet  this  sampling  rate  is  about  8 megahertz 
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The  2D  particle  display  system  is  a valuable  source  of 
data  to  LYC.  Data  is  collected  from  two  independent  systems 
(one  each  for  the  C130E  and  Lear jet)  and  recorded  on  magnetic 
tape  using  a nine- track  Pertec  recorder.  All  the  tapes  are 
processed  on  the  AFCRL  CDC  6600  computer. 

These  data  tapes  contain  two  types  of  records,  composed 
of  "fast"  and  "slow"  data.  The  fast  data  records  contain 
the  2D  particle  image  slices.  One  dimension,  the  columns, 
is  represented  by  the  32  diode  array,  the  other  dimension, 
the  rows,  by  time;  i.e.  one  row  represents  the  diode  status 
for  250  nanoseconds  of  time.  This  is  true  for  the  C130  air- 
craft only;  the  Learjet  row  of  data  occurs  every  125  nano- 
seconds because  of  the  faster  speed  of  the  jet. 

The  slow  data  records  occur  once  every  ten  seconds. 

These  records  contain  VCO  and  analog  information.  In  addi- 
tion, selected  ID  data  is  multiplexed  in  the  2D  buffer  and 
also  recorded  in  the  slow  data  records.  Note  that  the  ID 
data  appears  redundantly  (it  also  is  being  recorded  by  the 
Kennedy  seven  track  recorder) . 

The  exact  information  contained  in  these  records  is  dif- 
ferent for  each  aircraft.  All  the  programs  developed  for 
this  system  have  a variable  read-in  feature.  That  is,  each 
program  has  the  capability  of  reading  the  2D  tapes  recorded 
on  either  aircraft. 
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3.1  2D  PRE-PROCESSOR  (TWODEE) 

This  program  is  concerned  with  transforming  the  32  diode 
scans  into  actual  particles  or  crystals.  The  technique  de- 
veloped uses  the  real  time  principle  of  a one-pass  calculation; 
that  is,  while  the  particle  is  being  determined  certain  key 
parameters  are  also  being  calculated.  This  means  that  after 
the  initial  pattern  recognition  pass,  the  definition  of  the 
particle  is  fully  specified.  The  details  of  this  particle 
definition,  or  "signature",  are  explained  in  the  following 
sections . 

The  pattern  recognition  routine  employs  a string  tech- 
nique, where  a spring  is  defined  as  a series  of  consecutive 
occluded  diodes.  All  the  vital  information  required  for 
each  string  is  stored  in  one  computer  word,  including  (1) 
particle  number  or  identification,  (2)  scan  number,  (3)  be- 
ginning and  (4)  ending  diode  number,  and  (5)  a linkage  word. 

The  information  is  stored  as  five  12  bit  bytes  within  one 
60  bit  computer  word.  Routines  have  been  written  in  COMPASS 
to  insert  or  extract  a particular  byte  or  bytes  into  the 
computer  word.  Although  more  programming  is  required  using 
this  bit  manipulation;  it  is  more  desirable  than  using  5 
separate  words  requiring  5 times  as  much  storage. 

Each  new  string  found  is  checked  with  strings  from  the 
previous  scan.  If  the  string  is  adjacent  to  a previous  one 
(see  figure  7 ) it  is  considered  a string  of  the  same  parti- 
cle. The  new  string  identification  byte  must  be  copied  from 
the  previous  string  identification,  thus  identifying  the  new 
data  as  an  extension  of  the  previous.  The  old  string  linkage 
byte  will  have  to  be  set  to  point  to  the  new  string  word. 
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Two  strings  are  considered  adjacent  when  one  of  two  con- 
ditions exist:  (1)  If  two  strings  have  at  least  one  common 

diode  number  occluded  (the  left-hand  example  of  figure  7A 
shows  diode  #6  occluded  in  both  strings)  (2)  If  two  strings 
have  an  ending  occluded  diode  and  a beginning  occluded  diode 
within  one  number  of  each  other  (the  right-hand  example  of 
figure  7a  shows  the  beginning  diode  #3  in  the  first  string 
and  the  ending  diode  #2  in  the  second  string) . In  figure  7B 
neither  condition  exists,  and  the  strings  are  non  adjacent. 

2345678  123456789 

OOXXXOO  OOXXXXXOO 

OOOOXXO  XXOOOOOOO 

X = occluded  diode 
0 = lighted  diode 

Figure  7A:  Examples  of  non  adjacent  strings 

As  the  scanning  progresses  certain  particles  will  merge 
with  others  (see  the  stellar  example  in  figure  8 ) . The  rec- 
ognition scheme  using  this  definition  of  adjacent  strings 
accounts  for  this  perfectly,  and  the  particle  definition  in- 
formation will  be  combined.  After  all  1024  scans  have  been 
examined  the  particles  are  ready  for  the  next  processing 
step . 

OOOOOXXX  XXXOOOXXX 

XXXXOOOO  OOOOXOOOO 


Figure  7B:  Examples  of  non  adjacent  strings 
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For  those  4 strings  the  processor 
"sees”  3 different  particles 


3 particles  merge  at  this  time 

although  3 strings  are  present 
again,  they  are  all  adjacent  hence 
1 particle  remains 

particle  completed 
Particle  merging 

It  should  be  noted  that  only  completed  particles  are 
processed  at  this  time.  A completed  particle  is  defined  as 
a particle  completely  contained  within  the  1023  rows  in  the 
computer  memory.  If  a particle  has  a contribution  in  the 
1024th  scan  it  is  incomplete.  The  reasoning  for  this  incom- 
pletion is  that  the  particles  may  have  additional  strings  in 
the  next  record.  These  particles  are  saved  and  if  necessary 
merged  with  the  next  record.  In  either  case  they  become  the 
first  particles  processed  in  the  next  record . 


OOOOOOOOOOOOO 

looootfxtoo 

!OOQXXOOO 

:xqkxi6ooo 
oooo6^xxxxxx^6oooo 
oooooopKxxxCoooooo 

OOOOO^XXXXjMOOOOOO 

OOOC^XXO|Xxjb\<AOOOOO 

OOO^.^CKJQOO^XOOOO 

OOCKXOOOKXOOOKXpOO 

ooo^poooayoooayboo 

000000000000000000 


Figure  8 


A particle  is  said  to  be  defined  when  eight  fundamental 
parameters  about  it  are  known.  These  parameters  are:  area, 
perimeter,  volume,  horizontal  projection,  horizontal  Feret 
projection,  vertical  projection,  vertical  Feret  projection, 
and  longest  dimension.  Five  of  these  are  either  trivially 
calculated  or  a byproduct  of  the  string  technique  of  pattern 
recognition.  The  calculations  are  shown  in  the  following 
paragraphs . 


The  area  is  simply  the  total  number  of  diodes  occluded 


r 


by  a particle. 

AREA  = 1 number  of  diodes  occluded 

For  other  considerations  of  area,  see  Scientific  Report 
#1,  Determining  the  Volume  Represented  by  an  Irregularly 
Shaped  Cross-sectional  area,  1 Apr  75. 

The  perimeter  is  calculated  as  an  empirical  constant 
times  twice  the  sum  of  the  horizontal  projection  plus  the 
vertical  projection.  Figure  9 compares  the  actual  perim- 
eter with  the  calculated  perimeter  using  fb  = 1 . 


i 

i 
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where 


P = 2 (HP  + VP ) K j 

Kj  = a constant  for  particle  type  j 


Figure  9 : Perimeter  Calculation 

The  constant  fE  should  be  calculated  empirically  since  it 
may  be  too  high  for  soma  particle  types.  Until  that  time 


K.  should  have  a value  of  one 


Horizontal  Feret  projection  (hFP)  is  the  longest  singu 
lar  projection  in  the  horizontal  direction,  i.e.  the  direc- 
tion measured  by  the  diode  array  (along  the  row)  and  not  by 
time  (along  the  column) . This  horizontal  Feret  projection 
(figure  10)  is  calculated  by  taking  the  largest  occluded 
diode  number  less  the  smallest  occluded  diode  number  plus 
one. 


jf  = ending  diode  numbers 
jb  = beginning  diode  numbers 
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HFP  = 13  - 2 + 1 = 12  diodes 

Mathematically  the  vertical  Feret  projection  (VFP)  is 
similar  to  HFP,  the  difference  being  the  direction  of  mea- 
surement. Time  is  measured  by  scan  number  in  the  vertical 
direction.  The  VFP  calculation  is  the  last  scan  a particle 
appears  in  less  the  first  scan  a particle  appears  in  plus 
one . 


VFP  = MAX (nr)  - MIN (nr)  + 1 

where  nr  = scan  number 

Using  figure  10  the  VFP  is  calculated  as  ... 

VFP  = 12  - 2 + 1 = 11  diodes 

The  vertical  projection  (VP)  is  defined  as  the  sum  of 
the  partial  vertical  projections.  Figure  11  shows  this 
clearly.  The  equation  for  this  calculation  is 

VP  = PVPi  + PVP2 

where  PVP  = partial  vertical  projection 

The  individual  partial  projections  are  quite  cumbersome  to 
calculate . 

There  is  an  easier  way  to  calculate  VP  without  calcu- 
lating each  partial  vertical  projection.  VP  can  be  ascer- 
tained by  counting  the  the  number  of  strings  that  make  up  a 
particle.  This  can  be  seen  by  the  following. 


oxxxooooo 
ox  o oooo 
oxxxooooo 
OXXXOOOXO 
OCX  X V o v v o 

ocoxxxxxo 

ocoxxxxxo 

OOOXXXXOO 
OCOXXXXOO 
OOOO XX OOP 
OOOOOOOOO 


Figure  11:  Vertical  projection 


from  this  the  VP  is  calculated  by  adding  the  partial  projec 
tions : 


However,  this  is  precisely  the  number  of  strings  that  describe 
this  particle. 


The  horizontal  projection  (HP)  is,  of  course,  similar 
to  the  VP  except  for  direction.  A similar  analysis  can  be 
performed  to  calculate  HP.  The  particle  is  examined  for 
Strings,  not  in  the  "diode"  direction,  but  in  the  "scan" 
direction.  /fgain  the  number  of  strings  found  is  precisely 
the  horizontal  projection. 
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Volume  determination  is  an  area  where  we  have  done  ex- 
tensive research.  The  results  of  this  work  have  been  pub- 
lished in  a separate  scientific  report  entitled,  "Determina- 
tion of  the  Volume  Represented  by  an  Irregularly  Shaped 
Cross  Sectional  Area".  The  method  of  calculation  is  included 
in  this  report  but  the  reader  is  referred  to  the  original 
fora  complete  substantiation  of  the  conclusions.  The  method, 
simply  stated,  is  to  rotate  a two  dimensional  cross-sectional 
area  about  its  major  axis,  and  calculate  the  volume  generated. 
A correction  is  made  for  particle  irregularity  by  dividing 
the  area  into  two  halves,  one  above  and  one  below  the  major 
axis.  Actually  each  half  is  rotated  about  the  major  axis, 
and  a volume  is  calculated  for  each  half.  The  final  volume 
is  the  average  of  the  two  halves. 
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Consider  an  irregular  shape  in  2-space,  whose  volume  is 
desired : 


Firstly,  we  translate  the  axes  to  the  center  of  gravity  of 
the  geometric  shape.  Since  the  shape  is  formed  by  both 
interior  and  exterior  points,  all  points  are  considered. 
The  coordinates  of  the  center  of  gravity  (x„,y0)  are  given 
by 


xG  = l x./N 
k 

vG  - l ypn 

k 

where  x^,y^  are  the  points  along  and  inside  the  geometrical 
figure,  and  N is  the  number  of  x.,y^  points. 


We  then  form  the  new  noints  X . ,Y  . oiven  bv 
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Our  new  shape  is  unchanged  by  this  transformation.  The  cen- 
troid now  passes  through  the  origin. 


yt 


Figure  13:  Origin  at  centroid 

We  now  seek  the  major  axis  of  the  shape.  This  is  done  in 
order  to  rotate  about  the  axis  in  order  to  find  the  volume 
From  Vector  Mechanics  (see  Beer  & Johston,  STATICS  AND  DY- 
NAMICS, McGraw-Hill)  the  moments  of  inertia  for  the  shape 
are  given  by 


/ y2dA 
/ x2dA 
/ xydA 


in  the  discrete  case,  these  equations  become 


j 


- 


I - Z Y . 2 Am 
x 1 

I - Z X.2  Am 
Y i 

I - I X.Y.  4m 


I , I , I can  be  computed  at  the  same  time  , y are  done: 
x y xy  ij  (j 


Jy  = I (xi'XG)2  = Z xi2  - 2xGZxi  + Nxc 


Remembering  that 


1 Xi  = NXG 


= Z X.2  - 2Nx  2 + Nx  2 = Z x.2  - Nx  2 

1 Lj  (jj  1 Cj 


Similarly 


Ix  = 2 Yi2  - NyG2 
Xxy  = l (xi_XG)  (yi-yG} 

= 1 Xiyi  ~ XG  Eyi  " yG  Ixi  + NxGYC 
= Z xiYi  - NxGyG  - NxGyG  + NxGyG 

“ 1 Xiyi  ~ NxGyG 


thus  we  collect 


S = Z x. 
x 1 

S ? = Z x.  2 
x'-  1 


S = Z y. 

Y i 

S ? = Z v . 2 
T - i 


x-  y . 
x1 1 


*1 


(S  )2/N  I 
x ' x 

(Sx)  (Sy)/N 


(Sy) 2/N 


we  seek  an  angle  of  rotation  0 , such  that 

u.  = X.  cose  + Y.  sine 
11  1 

v^  = -X^  sin9  + cose 

would  force  the  new  uv  moment  to  vanish.  That  is 


I 

I 


uv 

uv 


l u . v . Am  = 0 

l l 

Am  l [X.  cose  + Y.  sine][-X.  sin9  + 

L l l l 


or 


l (Y.2  - X.2)  -5—  w2  + X.Y.  cos  20  = 0 

L l i 2 li 


substi tuting 


(1. 


- Iy) 


sin  26 


+ I 


xy 


cos  20  = 0 


or 


cose ]=0 


tan  20 
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After  rotating  through  the  angle  6,  we  have 


i 

I 

Figure  14:  Particle  rotated  to  major  axes 

Now  in  order  to  get  the  volume,  we  divide  the  area  into  two 
parts 

a)  all  points  (u^,v.)  which  lie  above  v = 0 

b)  all  points  (u^,v^)  which  lie  below  v = 0 

and  rotate  each  half  about  the  u axis.  The  resulting  volume 
will  be  the  average  of  the  two  rotated  areas. 

Pappus-Guldinus  second  theorem  states:  The  volume  of  a 
body  of  revolution  equals  the  generating  area  times  the  dis- 
tance traveled  by  the  centroid  of  the  area  while  the  body  is 
being  generated. 
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V = 2n  v A' 


where  v = the  centroid  component  in  the  v-direction 

A'  = area  being  rotated  (since  the  u-axis  bisects 
the  original  area,  this  area  is  one  half  ohe 
original  area;  A'  = A/2) 

each  centroid  v will  be  found  as 


vj  = E vi/k1 

kj  = number  of  points  v.  > 0 
v2  = E vi/k2 

k2  = number  of  points  < 0 


and  assuming  the  areas  Aj  and  A2 : 

Volume^  = 2tt  VjAj 
Volume2  = 2tt  v2A2 


or 


Volume  = 7T  (vjAj  + V2A2) 
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CONSIDERATION  OF  THE  AREA 

As  a result  of  the  considerations  investigated  in  Scien- 
tific Report  #1,  the  area  is  best  calculated  by  counting  oc- 
cluded diodes  above  and  below  the  major  axis.  When  an  occluded 
diode  center  is  conincindent  with  the  major  axis,  it  should 
be  counted  in  both  halves.  This  report,  and  the  scientific 
report  assumed  length  dimensions  of  1.  That  is,  the  distance 
between  diodes  multiplied  by  the  distance  between  successive 
diode  samplings  was  1.  The  final  Volume  should  be 

Volume  = tt  v N (Ap)  (Ax) 

where  Ap  = distance  between  diodes 

Ax  = sampling  distance,  determined  by  Aircraft 

speed  (m/sec)  multiplied  by  time  of  sampling 
rate  (sec) 

CONSIDERATION  OF  THE  ANGLE 


The  angle  of  rotation,  9,  is  given  by 

, 21 

9 = b tan  T--X^  j 

y x 

In  order  to  be  sure  that  this  rotation  will  yield  the  major 

axis  (rather  than  the  minor  axis)  the  arctangent  function 

will  be  required  to  produce  a result  between  -w  and  +ir  . If 

the  function  is  limited  to  a range  of  -j  and  j,  the  angle 

will  have  to  b°  modified  by  the  sign  of  I 

7 ^ xy 


pp 
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In  the  former  case,  the  resulting  arctangent  will  be 
between  -it  and  tt  ; Dividing  this  result  by  2 in  order  to 
obtain  9 will  yield 


IT 

2 


< 


9 


< 


IT 

2 


which  is  proper  rotation  to  insure  that  u will  be  the  major 
axis  and  v = u + v/2  the  minor  axis. 


Given  N points  (xi,yi)  i = 1,2,3,. ..,N 
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each  original  point  (x^y^) 


vi  = -c 2 (xi  - xQ)  + ci (yi  - 


E |vi  | 

count  all  v.  > 0 
1 - 

count  all  v.  < 0 
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The  longest  dimension,  LMAX,  will  be  found  to  lie  on 
the  major  axis.  The  volume  determination,  discussed  in  the 
previous  section,  calculates  this  axis  as  the  u-axis.  In 
addition,  all  the  x,y  particle  points  are  transformed  to 
u,v  coordinates.  The  longest  dimension  is  simply  the 
largest  u-value  less  the  smallest  u-value  plus  one.  Mathe 
matically  this  becomes 

LMAX  = MAX (u)  - MIN (u)  + 1 
Pictorially  this  is  represented  as 


Figure  15 : Longest  dimension 


3.2  2-D  PROCESSOR  ( KNOLL2D ) 

During  the  past  year  there  has  been  considerable  progress 
made  in  this  area.  A preliminary  version  of  program  KN0LL2D 
has  been  completed  and  is  currently  operational.  This  ver- 
sion categorizes  each  particle  by  size  class  along  two  sepa- 
rate dimensions,  maximum  length  (UMAX)  and  horizontal  Peret 
projection  (HFP) . Each  size  class  has  a length  of  two  diodes 
(i.e.  50  microns  for  the  cloud  probe  and  400  microns  for  the 
precip  probe) . The  number  of  observed  particles  in  each 
class  is  summed  over  1 second  intervals. 

Particles  that  touch  either  ending  diode  are  differen- 
tiated from  those  that  are  completely  enclosed  within  the 
32  diode  array.  This  requires  two  different  counts  to  be 
maintained;  one  for  the  total  number  of  particles,  and  one 
for  the  particles  completely  enclosed.  With  this  technique 
the  number  of  particles  that  would  have  been  ignored  had 
end-rejection  been  utilized  becomes  apparent.  A sample  of 
this  output  is  shown  on  the  following  page.  The  columns 
labeled  'W/0  ENDS'  are  the  counts  of  particles  NOT  touching 
an  ending  diode.  The  counts  under  the  column  labeled  'TOTAL' 
include  all  particles. 

An  additional  histogram  is  also  produced;  this  includes 
particles  categorized  by  area.  Both  the  total  and  non-re- 
jected  particles  are  shown  here  also.  The  sample  output  is 
included  on  the  following  pages. 
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Figure  16:  KN0LL2D  Sample  Output 
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Figure  17:  KN0LL2D  Sample  Output 
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The  following  tables  show  the  actual  channel  sizes, 
(center  diameter  and  upper  and  lower  limits)  and  cross  sec- 
tional area  (for  all  particles  and  the  non-rejected  particles 
only)  for  each  probe. 

CLOUD  PROBE 


:lass 

size (diodes ) 

lower 

upper 

center 

cross-sectional 

limit 

limit 

diameter 

end 

total 

(y) 

(y) 

(y) 

rej  ect 

(mm2 

) 

1 

1 



2 

12 . 5 



62.5 

37 . 5 

4.69 

5. 

0 

2 

3 

4 

62.5 

- 

112.5 

87 . 5 

17.50 

20. 

0 

3 

5 

6 

112.5 

- 

162.5 

137.5 

36.56 

45. 

0 

4 

7 

- 

8 

162.5 

- 

212.5 

187.5 

36.60 

48  . 

,8 

5 

9 

10 

212 . 5 

- 

262.5 

237.5 

33 . 55 

6 

11 

12 

262.5 

- 

312.5 

287.5 

30.50 

7 

13 

- 

14 

312.5 

- 

362.5 

337.5 

27.45 

[ 

8 

15 

- 

16 

362 . 5 

- 

412.5 

387.5 

24.40 

9 

17 

- 

18 

412.5 

- 

462.5 

437.5 

21.35 

10 

19 

- 

20 

462.5 

- 

512.5 

487.5 

18.30 

11 

21 

- 

22 

512.5 

- 

562.5 

537.5 

15.25 

12 

23 

- 

24 

562.5 

- 

612.5 

587.5 

12.20 

13 

25 

- 

26 

612.5 

- 

662.5 

637.5 

9.15 

14 

27 

- 

28 

662.5 

- 

"12.5 

687 . 5 

6.10 

15 

29 

- 

30 

712.5 

- 

762.5 

737.5 

3.05 

16 

31 

- 

32 

762.5 

— 

812.5 

787.5 

0.00 

i 

PRECIP  PROBE 

1 

1 

- 

2 

100 

- 

500 

300 

1536.0 

1638. 

4 

2 

3 

- 

4 

500 

- 

900 

700 

1495.2 

1708  . 

8 

3 

5 

- 

6 

900 

- 

1300 

1100 

1388.4 

t 

4 

7 

8 

1300 

- 

1700 

1500 

1281.6 

5 

9 

- 

10 

1700 

- 

2100 

1900 

1174.8 

6 

11 

- 

12 

2100 

- 

2500 

2300 

1068.0 

7 

13 

- 

14 

2500 

- 

2900 

2700 

961.2 

8 

15 

16 

2900 

- 

3300 

3100 

854.4 

9 

17 

18 

3300 

- 

3700 

3500 

747 . 6 

10 

19 

- 

20 

3700 

- 

4100 

3900 

640.8 

11 

21 

- 

22 

4100 

- 

4500 

4300 

534.0 

12 

23 

24 

4500 

- 

4900 

4700 

427.2 

13 

25 

26 

4900 

- 

5300 

5100 

320.4 

14 

27 

28 

5300 

- 

5700 

5500 

213.6 

15 

29 

- 

30 

5700 

- 

6100 

5900 

106.8 

16 

31 

- 

32 

6100 

- 

6500 

6300 

0.0 

1 

: 
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The  remaining  analysis  done  in  this  program  uses  the 
particle  counts  with  end-rejection.  The  reader  should  keep 
in  mind  that,  unlike  the  ID  system,  the  number  of  particles 
ignored  by  this  feature  is  readily  available,  by  comparing 
the  tables  previously  discussed. 


Two  new  parameters  have  been  defined  and  are  calculated 
for  each  particle;  these  are  1)  average  projection  ratio  and 
2)  equivalent  circle  ratio.  The  average  projection  ratio  is 
a measure  of  a particle's  contour.  That  is,  the  more  closely 
'stellar  shaped'  a particle  is,  the  higher  its  average  pro- 
jection ratio.  Mathematically  the  formula  is 


APR 


1 

2 


(J* 

lHFP 


VP 

VFP 


) 


where 

and 


HP  and  VP  are  total  projections 
HFP  and  VFP  are  Feret  projections 


The  equivalent  circle  ratio  is  a measure  of  a parti- 
cle's deviation  from  a circle.  The  equation  for  equiva- 
lent circle  ratio  (ECR  = P/T 4ttA)  is  derived  as  follows: 

ECR  _ P/A  where  the  c subscript  indicates 

P /A  equivalent  circle  perimeter 

and  area 


but  A = A 

c 

P 

ECR  = ^ 
c 

P = 
c 


where 


itD 


These  two  parameters  are  categorized  into  sixteen  size 
classes  and  the  number  of  particles  in  each  class  is  shown. 
The  size  classes  for  the  average  projection  and  equivalent 
circle  ratios  are  shown  in  the  following  tables. 


class 

apr  : 

Limits 

ECR  : 

Limits 

1 

1.00 

1.10 

1.00 

1.20 

2 

1.10 

1.20 

1.20 

1.40 

3 

1.20 

1.30 

1.40 

1.60 

4 

1.30 

1.40 

1.60 

1.80 

* 

5 

1.40 

1.50 

1.80 

2.00 

\ 

6 

1.50 

1.60 

2.00 

2.20 

7 

1.60 

1.70 

2.20 

2.40 

8 

1.70 

1.80 

2.40 

2.60 

9 

1.80 

1.90 

2.60 

2.80 

10 

1.90 

2.00 

2.80 

3.00 

11 

2.00 

2.10 

3.00 

3.20 

12 

2.10 

2.20 

3.20 

3.40 

13 

2.20 

2.30 

3.40 

3.60 

14 

2.30 

2.40 

3.60 

3.80 

15 

2.40 

2.50 

3.80 

4.00 

16 

2.50 

2.60 

4.00 

4.20 

65 


Four  sets  of  matrices  are  utilized  to  further  describe 
the  particle  distribution.  The  matrices  display  the  particle 
counts  by  class  size  in  the  following  manner: 

1)  UMAX  and  HFP 

2)  UMAX  and  AREA 

3)  UMAX  and  APR 

4)  UMAX  and  ECR 

Each  matrix  is  calculated  twice;  once  for  the  cloud  probe 
and  once  for  the  precip  probe.  Examples  of  these  matrices 
are  shown  on  the  following  four  pages. 

Another  area  where  development  has  begun  is  in  automa- 
tic particle  typing.  A display  of  the  average  projection 
and  equivalent  circle  ratios  in  a common  matrix  provides 
insight  into  the  composition  of  the  cloud  particles.  This 
matrix  is  divided  into  sections  of  known  particle  types. 

The  relative  percentages  of  each  particle  type  present  are 
easily  calculated.  Examples  of  both  the  theoretical  and 
actual  matrices  are  shown  in  the  following  pages.  In  the 
actual  case,  the  data  collected  was  in  predominantly  rain, 
as  borne  out  by  the  matrix  display. 

The  programming  effort  described  in  this  section  is  by 
no  means  complete.  The  work  that  has  been  detailed  here  is 
really  just  a beginning  of  the  development  required  to  glean 
the  maximum  amount  of  information  from  this  data  source.  It 
is  a project  that  will  require  continual  modification  and 
communication  between  DPSI  and  the  LYC  scientists. 
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Figure  19:  KN0LL2D  Sample  Output 


Figure  21:  KN0LL2D  Sample  Output 


EQUIVALENT  CIRCLE  RATIO 


APR  AND  lITUr/4L£vll  ClRCLE  iAIIU  CLJJO  P* 


Figure  23:  KNOLL2D  Sample  Output 
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3.3  2D  UTILITY  PROGRAM  (KN2UTIL) 

Because  of  the  voluminous  amount  of  data  collected  by 
the  2D  particle  display  system,  a few  considerations  should 
be  made  before  any  processing  is  done.  For  example,  only 
data  between  the  specified  pass  times  should  be  processed. 
Also  any  data  that  appears  faulty  or  even  marginal  should 
be  eliminated.  This  program  provides  the  verification  re- 
quired for  these  determinations.  The  program  has  been  writ- 
ten to  employ  two  different  output  media.  Each  output  de- 
vice has  its  own  advantages  and  a definite  place  in  the  pro- 
cessing stream. 

The  6600  line  printer  output  is  one  means  of  displaying 
the  observed  particles.  This  printer  should  be  utilized 
when  there  is  a limited  number  of  records  to  be  examined, 
or  when  fast  inspection  of  the  data  is  vital  to  LYC . One 
might  expect  to  get  two  or  three  of  these  computer  runs  per 
day  using  the  line  printer. 

This  limitation  is  due  to  the  long  record  length,  and 
it  becomes  impractical  to  print  many  records.  Each  record 
contains  1024  scans  of  the  32  diode  array,  yielding  approxi- 
mately 20  computer  pages  per  record.  A maximization  scheme 
was  developed  that  allows  three  records  to  be  listed  using 
20  pages.  Even  with  this  technique  to  reduce  output  lines, 
a 30  record  listing  produces  200  pages.  Therefore,  although 
the  line  printer  output  offers  faster  results,  its  use  is 
limited. 

Another  method  of  displaying  this  data  is  utilizing  the 
35  mm  film  plotter  available  at  the  computer  center.  There 
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are,  again,  advantages  and  disadvantages  of  this  technique. 
The  primary  advantage  of  using  this  medium  is  that  one  full 
record  may  be  displayed  on  a single  35  mm  frame.  Utilization 
of  the  film  copier  availabe  at  LYC  provides  any  hard  copies 
desired.  In  addition  storage  of  the  film  rolls  is  much  more 
desirable  than  storing  cumbersome  computer  listings.  The 
only  drawback  to  this  method  is  turn-around  time.  The  com- 
puter center  requires  one  full  day  for  film  processing. 

To  summarize,  the  following  results  are  produced  by 
this  program: 


A tape  summary  of  all  data  recorded  that  includes 

a.  record  type  (slow  or  fast) 

b.  record  length 

c.  record  number  (absolute  and  by  type) 

d.  record  time 

A data  listing  by  record  type  as  specified 

a.  slow  data  - line  printer 

b.  fast  data  (on  selected  device) 

1.  line  printer 

2 . 3 5 run  f i lm 


Because  the  two  record  types  are  interleaved,  each  re- 
cord to  be  listed  may  be  specified  in  one  of  two  ways,  either 
by  overall  record  number  or  by  sequence  number  within  type 
(fast  or  slow) . If  a tape  consists  of  the  following  records, 
the  record  count  is  as  shown. 


S S S F F S FI  F S S F F F F S | S F F S F 

Record  number  1 2 3 4 5 6 7 8 9 10  11  12  13  14  15  16  17  18  19  20 

Slow  sequence  #123  4 56  78  9 

Fast  sequence  # 1234  5678  9 10  11 
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Case  A. 


To  dump  all  of  the  above  records  either  of  the 
following  methods  should  be  used. 


1.  records 


slow  records 


fast  records 


1-20 


1-11 


Case  B. 


To  dump  the  records  indicated  by  the  above  eec- 
tangle  use  either: 


1.  records 


8-15 


2.  slow  records  5-7 

and 

fast  records  4-8 

Records  to  be  dumped  are  input  via  a $SS  card  in  stan- 
dard namelist  format.  The  record  number  technique  uses  pre- 
fix A;  the  slow  record,  prefix  S;  and  the  fast  record  prefix 
F.  The  control  variable  for  the  $SS  card  is  REC  prefixed 
with  this  character,  and  suffixed  with  a B or  E indicating 
beginning  or  end  of  dump.  For  the  previous  two  cases  the 
$SS  card  should  be: 


Case  A. 


@$SS  ARECB=1 , ARECE=2  0 $END 


@$SS  SRECB=1 , SRECE=9 , FRECB=1 , FRECE  = g $END 


Case  B. 


@$SS  ARECB=8,ARECE=15  $END 


where 


@$SS  SRECB=  5 , SRECE=7 , FRECB=4 , FRECE=  6 $END 
@ = blank  column 


LEARJET  DATA  COLLECTION  SYSTEMS 


LYC  receives  data  tapes  from  two  independent  systems 
on  board  the  Learjet  36C  high  altitude  aircraft.  These  are: 
the  PMS  ID  particle  sizing  system  and  the  PMS  2D  particle 
display  system.  The  2D  system  is  identical  to  that  on  the 
C130E  aircraft,  with  the  exception  of  a few  minor  VCO  loca- 
tions. As  discussed  in  the  previous  section  the  current  2D 
programs  will  read  and  process  this  data  without  any  diffi- 
culty. The  PMS  ID  system  did  require  an  additional  program- 
ming effort. 

There  are  two  types  of  ID  data  received  at  LYC;  processed 
and  unprocessed  (raw) . The  raw  data  tapes  have  two  major  dif- 
ferences from  the  tapes  produced  on  the  C130  aircraft.  The 
PMS  ID  tape  produced  on  the  Learjet  is  9 track,  800  bpi. 

Data  is  collected  at  the  rate  of  64  words  and  blocked  with 
32  seconds  per  record.  (The  tapes  produced  on  the  C130  air- 
craft also  contain  64  words  per  second;  however  they  are 
blocked  with  4 seconds  per  record.)  The  word  structure  on 
the  Learjet  tapes  consist  of  four  4-bit  BCD  characters,  for 
a word  length  of  16  bits.  (The  C130  word  structure  consists 
of  four  6-bit  characters,  for  a word  length  of  24  bits.)  Also 
the  recording  system  utilized  on  the  Learjet  is  a nine-track 
Kennedy  1600/360  recorder,  (the  C130  system  uses  a Kennedy 
seven-track  tape) . 

With  the  addition  of  two  new  programs  and  modification 
of  a third,  the  complete  Learjet  data  collection  system  has 
been  incorporated  into  the  LYC  data  processing  stream. 
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4.1  ID  REFORMATTER  (LEARPMS) 

The  processing  of  data  is  a rather  complete  system  when 
the  origin  of  the  data  is  the  C130  aircraft.  Learjet  data 
is  similar  to  and  should  be  analyzed  in  the  same  way  as  the 
C130  data.  It  is  most  advantageous,  then,  to  convert  the 
raw  Learjet  data  to  a form  readable  by  the  present  ID  pro- 
cessing system  so  that  the  data  can  be  analyzed  as  if  it 
were  from  the  C130.  This  decision  has  two  advantages.  Firstly, 
it  eliminates  many  modifications  to  KN0LL1D  which  would  only 
serve  to  further  complicate  this  program.  Secondly  this  de- 
cision provides  for  a backup  tape  to  be  produced.  It  should 
be  noted  that  there  is  no  other  backup  tape  available  from 
the  Learjet  aircraft. 

Program  LEARPMS  has  basically  two  functions.  It  will 
read  a nine  track  32  second/record  tape,  perform  the  neces- 
sary reblocking  and  produce  a seven  track  4 second/record 
type.  When  this  new  tape  is  created  it  will  be  input  di- 
rectly to  the  existing  KNOLL ID  program. 

This  program  has  been  written  and  completely  checked 
out.  The  program  requires  no  input  deck;  it  simply  reads 
the  input  tape  and  produces  a reformatted  output  tape, 
suitable  for  KNOLLlD  processing. 
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4.2  ID  POST  PROCESSOR  (HIAC1D) 

The  second  source  of  Lear jet  ID  data  LYC  receives  is 
the  contractor-processed  data.  This  is  recorded  on  seven- 
track  magnetic  tape  in  IBM  floating  point  format. 

Program  HIAC1D  was  written  to  read  this  tape;  reformat 
the  data  to  CDC  floating  point;  perform  basic  editing,  if 
desired;  calculate  the  normalized  particle  number  density; 
average  the  results  over  any  given  interval;  calculate  the 
median  volume  diameter  for  each  probe;  list  the  results  in 
an  output  similar  to  KNOLL ID ; and  produce  an  output  tape 
suitable  for  processing  with  program  RAPP. 

This  output  tape  is  identical  to  the  RAPP  tape  produced 
by  program  KN0LL1D.  It  allows  the  Lear jet  data  to  be  input 
to  the  Radar  and  Aircraft  Plotting  Program  (RAPP  see  sec.  5) 
and  correlated  with  the  radar  data  received  at  LYC. 
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4.3  ID  PLOTTER  (KNPLT1D) 

An  option  has  been  written  for  the  standard  PMS  ID 
plotting  program  (KNPLT1D)  that  allows  this  data  to  be 
plotted.  When  specified,  this  program  reads  the  Learjet 
tape;  performs  the  necessary  conversions  and  calculations; 
and  produces  a file  similar  to  the  plot  tape  created  in 
program  KNOLLlD.  The  program  then  produces  a set  of  plots 
similar  to  those  produced  from  the  data  collected  by  the 
C130  aircraft. 


RADAR  DATA  ANALYSIS 


During  each  C130  data  collection  mission,  the  aircraft 
is  tracked  by  a ground  based  weather  radar  facility.  This 
radar  samples  the  environment  in  a known  volume  of  space 
just  before  the  aircraft  passes  through  it.  It  provides  a 
measure  of  the  radar  reflectivity  factors  for  the  ice  cry- 
stals and  water  droplets  being  observed  by  the  aircraft. 

The  processing  stream  involved  to  utilize  this  data  re- 
quires two  programs  to  be  run.  The  first  program,  SPANDAR, 
is  a data  display  and  reformatting  program.  It  also  produces 
an  output  tape  to  be  used  with  the  analysis  program,  RAPP. 
Program  RAPP  performs  the  analysis  and  correlations  between 
the  radar  and  aircraft  data.  Details  of  these  two  programs 
are  shown  in  the  following  sections. 


RADAR  DATA  DISPLAY  (SPANDAR) 


The  radar  link  tapes,  received  at  LYC,  contain  the 
meteorological  information  as  seen  by  the  C130  aircraft 
during  its  flight.  Unfortunately  these  tapes  are  produced 
at  the  Applied  Physics  Laboratory  in  Baltimore,  where  an 
IBM  installation  is  employed,  and  the  data  is  in  IBM  float- 
ing point  format.  A program  has  been  developed  to  pre- 
process  these  tapes.  That  is,  a routine  was  written  to 
convert  the  format  from  32  bit  IBM  360  floating  point  to 
60  bit  CDC  6600  floating  point.  The  routine  also  converts 
8 bit  EBCDIC  character  data  to  6600  display  code. 

During  the  past  year  the  blocking  algorithm  utilized 
by  APL  was  changed  unexpectedly.  To  account  for  this,  pro- 
gram SPANDAR  had  to  be  changed;  the  program  now  accepts  any 
blocking  format  for  a single  tape.  However,  an  automatic 
dump  procedure  has  been  incorporated  into  the  program.  If 
the  blocking  format  keeps  changing;  the  auto-dump  mode  is 
selected;  the  block  in  error  is  dumped  hexadecimally  and  the 
program  terminated.  This  facilitates  any  future  changes  as 
a result  of  different  blocking  techniques. 

SPANDAR  produces  a tape  summary  that  includes,  for  each 
header  record,  the  following  information:  aircraft,  pass, 
date,  height,  and  the  data  start  and  stop  times.  An  output 
tape  is  also  produced.  This  tape  is  used  as  the  input  radar 
tape  for  program  RAPP.  These  two  functions  are  produced 
with  every  SPANDAR  run. 

The  program  has  three  additional  options,  selected  by 
control  cards  at  the  users  request.  These  include:  1)  a line 
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printer  plot  of  radar  ref lectivitiy  (in  DB)  as  a function 
of  time  and  radar  slant  range,  2)  a listing  of  the  input 
tape,  and  3)  a punched  deck. 


t 


5.2  RADAR  AND  AIRCRAFT  DATA  COMPARISON  (RAPP) 

Program  RAPP  is  the  means  of  comparing  the  actual  liquid 
water  content,  as  derived  by  the  particle  measurements  ob- 
tained on  the  C130  aircraft,  against  the  total  mass  of  the 
observed  particles,  as  inferred  from  the  measured  radar  re- 
flectivity. The  program  uses  two  input  tapes,  one  containing 
aircraft  data,  the  other  containing  radar  data.  The  aircraft 
generated  data  is  produced  as  an  output  tape  from  program 
KNOLLlD.  The  radar  tape  produced  is  a result  of  program 
SPANDAR.  The  variables  of  concern  are  time,  liquid  water 
content  (mass),  and  radar  reflectivity.  Liquid  water  con- 
tent is  itself  a function  of  the  selected  probe;  the  program 
will  use  the  liquid  water  content  for  either 

(1)  scatter  probe 

(2)  cloud  probe 

(3)  precip  probe 

(4)  total  of  cloud  plus  precip  probes 

A full  set  of  plots  (section  5.2.5)  is  produced  for  each 
pass.  In  addition,  the  radar  and  aircraft  data  is  listed 
along  with  the  correlation  and  curve  fitting  coefficients. 

5.2.1  DATA  RECONSTRUCTION 

KNOLLlD  outputs  its  aircraft  data  every  n seconds  (where 
n is  a parameter  of  the  program) , while  SPANDAR  outputs  data 
every  second.  In  order  to  make  a better  comparison  of  the 
data,  the  aircraft  data  must  be  known  for  every  second.  This 
KNOLLlD  output  is  the  arithmetic  mean  of  the  n-point  smoothed 
data,  and  reported  at  the  beginning  of  every  n-second  time 
interval . 

In  order  to  reconstruct  one  second  data,  two  methods 
come  to  mind  (let  us  assume,  that  for  the  sake  of  this  dis- 
cussion , n = 4 ) : 


fed : 


(a)  repeat  the  same  data  for  each  second  four  times, 
and,  by  taking  3 such  sets,  find  the  least  square 
quadratic  to  the  12  points;  evaluate  the  quadra- 
tic at  each  second. 

(b)  find  the  quadratic  exactly  fitting  3 points,  each 
of  which  is  4 seconds  apart;  evaluate  the  quadra- 
tic at  each  second. 

The  first  method  involves  generation  of  Ex,  Ex2,  Ex3, 
Ex4,  Ey,  Exy,  Ex2y  for  the  equally  spaced  data.  One  could 
develop  the  matrix 


n 

EX 

Ex2 

Ex 

Ex2 

Ex3 

Ex2 

Ex3 

EX4 

and  get  its  inverse. 

since 

the  x-values  never  change. 

Then 

multiplying  this  inverse  by  the  vector 

Exy 

lx2y 

would  yield  the  results  c,  b,  a in 
y = ax2  + bx  + c 

The  second  method  is  more  easily  solved,  using  Newton's 
Forward  formula  over  the  three  equally  spaced  points: 

Y = YO  + u (y i -y o ) + 


(y2-2y!+  y0) 


where 


u 


< 


The  time  variable,  x,  will  take  on  the  values,  in  seconds, 
from  the  initial  time  period, 

T = .5,  1.5,  2.5,  3.5, ... , 11.5 

where  the  values  of  y are  known  at 
T = 2,  6,  10 

Substituting  these  values  in  the  equation,  we  get 


T-2 


y = yo+±Ti(yi-yo)  + (±4^)  (^p) 


-2.  ,t-6. ,y2-2yi+y° 


evaluating  at 


T = .5,  1.5 


we  get 

y.5  = y0~ • 375  (yx -y0 )+ (-.  375)  J)  (y2~2y1+y0) 

y>5  = y o - • 375 (y ! -y0 )+. 2578125 (y2~2y1 +y0 ) 
y>5  = 1.6328125y0-. 89062 5yi+. 257812 5y 2 


and 


*1.5  - y0--125(y,-y0H-<-. 125X11415,  (Y2.2yi+yo) 
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T-TT-  17’ 


: 


Y1.5  = Yo~ • 125 (y1-y0)+. 0703125 (y2-2y1+y0) 
y1  5 = 1. 195312 5y0-. 26562 5yi+. 070312 5y 2 


From  the  first  two  points  y(.5),  y { 1 . 5 ) , we  can  generate 
the  entire  set  of  y by  a table  of  differences,  namely 


1.5 

2.5 


. 5 


1.5 


2 . 5 


Ay 


. 5 


A 2y 


. 5 


Ay 


1.5 


The  final  column  A:y  must  be  constant  since  a quadratic 
fits  all  the  data,  so  each  subsequent  y can  be  calculated 


yk+ . 5 


A 2y 


iyk-1.5  + yk- 


5 


let 


A2y  5 = A2,  a constant 


yk+ . 5 = A2  + (yk-.5"yk-1.5)  + yk- . 5 


'k+.5  L‘  + 2yk- . 5 yk-l . 5 


for 

k = 2,  3,  4,  5, . . . , 12 

A'"  can  be  verified  to  be  (using  T = 2.5) 


(yn-2y,+y,) 
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Procedure  for  interpolation: 


1.  calculate  A2  = . 0625 (70-27!+72) * 


2.  calculate  yx  and  y2i 


7i  = 1.63281257Q-. 8906257J+. 257812572 


72  = 1.19531257Q-. 2656257!+. 070312572 


3.  set  k = 3 


4 .  evaluate 


Yk  = 2yk-l  ’ yk-2  + A2 


5.  increase  k b7  1;  if  k > 12  stop  otherwise  return  to  step 


'note : 


.0625 

1.6328125  = 
.890625  = 

.2578125  = 
1.1953125  = 
.265625  = 

.0703125  = 


1/16 

209/128 

114/128 

33/128 

153/128 

34/128 

9/128 


This  procedure  must  be  compared  with  the  least  squares 
algorithm  in  order  to  determine  the  more  desirable  fitting 
technique.  Three  examples  were  performed,  and  referred  to 
as  examples  1,  2 and  3. 

Example  1:  at  time  =2,  y = 4 

at  time  = 6,  y = 6 
at  time  = 10,  y = 9 

the  computed  results  follow: 


T 

Nominal 

Newton's  method 

Least  squares 

. 5 

4 

3 .5078125 

3.549454 

1.5 

4 

3.8203125 

3.949051 

2.5 

4 

4.1953125 

4.380619 

3.5 

4 

4 . 6328125 

4.844156 

4.5 

6 

5.1328125 

5. 339660 

5.5 

6 

5 . 6953125 

5.867132 

6.5 

6 

6.3203125 

6.426573 

7.5 

6 

7 . 0078125 

7 . 017982 

8.5 

9 

7.7578125 

7.641358 

9.5 

9 

8 . 5703125 

8.296703 

0.5 

9 

9.4453125 

8 . 984016 

1.5 

9 

10.3828125 

9.703297 

Figure 

24  is  the 

comparison  of  these 

plots  for  example  1. 

The 

dashed  line  on 

the  plot  is  for  Newton's  method?  the  solid 

for 

least 

squares . 

The  + indicates  the 

given,  or  nominal 

points . 

Correspondingly  the  first  equation  is  for  the  Newton 

89 


l 


I 

W 


| . 


I 


tecnnique,  the  second  equation  for  least  squares.  Note  that 
the  dashed  line  passes  through  the  midpoint  of  the  horizontal 
line  connecting  each  set  of  four  points.  It  appears  that 
either  curve  would  be  satisfactory  in  this  case. 

Example  2:  at  time  = 2,  y = 4 

at  time  = 6,  y = 8 

at  time  = 10,  y = 12  (straight  line) 


T 

Nominal 

Newton's  method 

Least  squares 

. 5 

4 

2.5 

3 . 076923 

1.5 

4 

3.5 

3.972028 

2.5 

4 

4 . 5 

4.867133 

3.5 

4 

5.5 

5.762238 

4.5 

8 

6.5 

6.657343 

5.5 

8 

7.5 

7.552448 

6.5 

8 

8.5 

8 . 447552 

7.5 

8 

9.5 

9.342657 

8.5 

12 

10.5 

10.237760 

9.5 

12 

11.5 

11.132870 

10.5 

12 

12.5 

12.027970 

11.5 

12 

13.5 

12.923080 

Again, 

the  dashed 

line  on 

the  plot  is 

for  Newton's  method 

and  the 

solid  is 

for  the 

least  squares 

The  curves  are 

again  relatively  close,  and  there  appears  to  be  little  mo- 
tivation for  selecting  one  method  in  preference  to  the  other. 

Figure  25  compares  these  two  plots. 
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Example  3: 


Nominal 


10.5 

11.5 


at  time  =2,  y = 4 
at  time  =6,  y = 9 
at  time  = 10,  y = 6 

Nev/ton's  method 


. 0625 

2.8125 

5.0625 

6.8125 

8.0625 

8.8125 

9.0625 

8 .8125 

8 . 0625 

6.8125 

5. 0625 

2.8125 


Least  squares 


2.758242 
4 . 260740 
5.507492 
6.498501 
7.233766 
7.713286 
7. 937062 
7.905094 
7 . 617382 
7.073926 
6.274725 
5.219781 


See  figure  26  for  the  comparison  of  these  curves.  Here 
we  see  an  interesting  phenomenon.  The  dashed  line  (Newton's) 
fits  the  data  considerably  better  than  the  least  square  curve. 
Remembering  that  these  given  points  have  already  been  smoothed 
by  averaging,  it  appears  that  the  least  square  technique  is 
smoothing  even  further.  That  is,  the  least  square  technique 
is  smoothing  out  the  peaks  and  valleys,  eliminating  the  ori- 
ginal character  of  the  data. 

This  example  leads  to  the  conclusion  that  Newton's 
method  is  the  more  desirable  tecnique.  It  is  noteworthy 
that  it  is  also  easier  and  thus  more  efficient. 


Although  this  technique  works  well,  a recent  flight  (27  MAR 
76)  did  yield  erroneous  data.  Because  the  LWC  and  Z averages 
were  so  low,  some  calculated  points  along  the  curve  entered 
the  negative  domain.  This  was,  of  course,  unacceptable  and 
a new  technique  was  considered. 

Newton's  Forward  formula,  however,  did  work  well  with 
good  data.  The  same  technique  was  attempted  using  the  logs 
of  each  point,  rather  than  actual  data.  After  the  interpo- 
lation, antilogs  were  taken,  and  the  results  were,  of  course, 
always  positive.  This  modification  worked  well  and  is  currently 
being  used.  The  steps  used  are  shown  below: 

For  three  given  points  (tj,yi),  (t2,y2),  ( 1 3 r y 3 ) and  an 

average  interval  m 

1.  take  logs 

Yj  = LOG10(y1) 

Y2  = LOG1Q(y2) 

Y 3 = LOG1Q(y3) 

2.  calculate  Y for  each  time  t,  from  t = ti-m/2  to 
t = t3+m/2  using  . . . 

Yt  - Yj  + (t-j)[Y2  - Yj  + <t~)  (Y3  -2Y2  + Yj)|]/m 


3.  take  antilogs 
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5.2.2  AIRCRAFT-RADAR  CORRELATION 

The  radar  data,  as  has  been  explained,  is  the  result  of 
locking  the  signal  at  a fixed  distance  in  front  of  the  air- 
craft. This  means  that  the  radar  time  for  the  particular 
data  is  early,  and  the  time  has  to  be  changed.  By  examining 
the  radar  data  of  azimuth,  elevation  and  range,  one  can  de- 
termine the  ground  speed  over  a desired  time  interval.  Then 
one  can  calculate  the  time  lag,  t 

t = d/s 


where  d is  the  "lock-on-distance"  and  s is  the  ground  speed. 
This  value  of  t is  added  to  the  reported  value  of  t for  the 
radar  data.  However,  because  of  the  nature  of  the  radar  in- 
strumentation, there  might  still  be  a timing  error  of  a few 
seconds.  In  order  to  minimize  this  effect,  a cross-correla- 
tion is  performed  between  the  aircraft  z and  the  radar  z. 

This  correlation  is  done  with  At  = -3,  -2,  -1,  0,  +1,  +2,  +3 
and  the  At  corresponding  to  the  maximum  correlation  is  chosen 
as  an  additional  time  offset. 

The  equations  used  (REF:  INTRODUCTION  TO  MATHEMATICAL 
STATISTICS,  Hoel,  pp  117-122)  is  shown  below. 


n 

l [ (x  ■ --x)  (y  . -y)  ] where: 


r = 


i=l 


n S £ 
x y 


S = standard  deviation 
r = correlation  coefficient 


r = 


Exy  - nxy 


/[  lx2~nx2  ] [ Ey2-ny2  ] 


J 


/[n£x2- (Ex) 2][nEy2-(Ey)2] 

LEAST  SQUARE  CURVE  FITTING 

A least  squares  fit  to  the  data  is  made  for  plots  1, 

2 and  5 from  list  on  page  96.  Since  the  choice  of  variables 
for  x and  y is  rather  arbitrary,  two  least  square  lines  are 
calculated: 


y = ax  + b 
x = cy  + d 

and  the  RMS  (root  mean  square)  error  is  calculated  and  re- 
ported for  each.  The  fit  with  the  lesser  RMS  error  is  chosen 
as  the  fitting  function.  The  formulas  for  the  fitting  func- 
tions and  RMS  are: 


x = £x/m 

y = Ey/m 


(m  points) 


a = 


b = 


Exy  - yEx 
Ex2  - xEx 

y - ax 


c = Exy  - xEy 
£y2  - yEy 

d = x - cy 


r = Ey2  + a2 Ex2  - 2aExy 
1 m 


b2 


_ Ex2  + c2Ey2  - 2c Exy  _ ,2 


m 


1 


At  the  request  of  Mr.  Vernon  Plank  a filtering  routine 
was  added  to  program  RAPP.  The  routine  allows  selected  radar 
data  (Z)  to  be  omitted  if  it  falls  below  the  minimum  detect- 
able radar  Z.  This  minimum  detectable  signal  is  not  a con- 
stant but  a function  of  radar  slant  range  and  a calibration 
constant.  The  equation  used  is  as  follows: 


Z . (DBZ ) = DB,  + 20  LOG. n (r ) + corr 
mm  1 10 

where 

DB!  = calibration  constant  at  1 km 
r = radar  slant  range  (km) 
corr  = particle  type  correction  6.5 
for  ice,  0.0  for  water 

This  function  is  also  plotted  on  the  radar  Z histogram  to 
graphically  depict  the  data  to  be  ignored. 

5.2.5  SAMPLE  PLOTS 

The  list  of  plots  below  are  included  on  the  following 
pages. 

1.  A/C  MK  (1000*M/vTz)  vs.  RADAR  Z (Log-Log) 

2.  A/C  M vs.  A/C  Z (Log-Log) 

3.  RADAR  and  A/C  Z vs.  TIME  (Semi  log) 

4.  A/C  MK  and  RADAR  Z vs.  TIME  (Semi  log) 

5.  A/C  Z vs.  RADAR  Z (Log-Log) 

6.  FLIGHT  TRACK  (FOR  KWAJALEIN  DATA  ONLY) 
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6.  AIRBORNE  DATA  COLLECTION  SYSTEM 


The  airborne  data  collection  system  utilized  by  LYC  is 
composed  of  a Digital  Equipment  Corporation  PDP  8/E  central 
processor  with  an  8K  core  memory.  The  following  peripheral 
devices  are  also  included: 

1.  twin  dectapes 

2.  Dec  magnetic  tape  unit 

3.  GE  high  speed  printer/keyboard 

4 . Tektronix  CRT/keyboard 

5.  PMS  ID  data  buffer 

6 . A to  D converter 

Figure  27  on  the  next  page  shows  a block  diagram  of  the 
complete  data  collection  systems  on  board  the  C130A  and  E air- 
craft. 

Figure  28  on  the  next  page  following  depicts  the  PDP  8/E 
and  its  associated  peripherals  on  board  the  C130E. 

There  also  exists  a PDP  8/E  at  LYC  which  is  similar  to 
the  airborne  PDP  8/E,  and  is  used  primarily  for  program  de- 
velopment. The  configuration  of  peripheral  hardware,  however, 
on  the  LYC  computer  is  not  the  same  as  the  airborne  machine. 

At  LYC,  there  is  neither  an  A to  D converter  nor  a GE  I/O 
High  speed  printer/keyboard.  The  DEC  writer  is  used  for 
inputting  numbers,  and  for  the  printing  of  the  results.  Thus 
no  VCO  or  hard-wired  input  comes  into  the  computer  directly. 

At  LYC,  inputs  which  are  necessary  for  debugging  are  simulated 
by  patch  boards,  or  the  reading  of  output  tapes  back  into  the 
computer . 


Figure  27:  Data  flow  on  the  C130A  and  C130E  aircraft 


PMS  ID 
BUFFER 


Figure  28:  Airborne  PDP  8/E  with  Associated  Peripherals 
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In  addition,  there  are  significant  differences  in  the 
IOT  (input-output  transfer)  commands.  The  IOT's  are  the 
commands  the  programmer  uses  in  order  to  obtain,  make  ready 
or  output  parameters  from  or  to  core  locations.  The  reader 
should  note  that  certain  IOT's  are  not  available  on  the  LYC 
computer  (PMS  and  VCO  inputs) . More  importantly,  the  pro- 
grammer is  cautioned  that  a particular  IOT  will  execute  quite 
differently  on  the  two  computers.  Also,  both  computers  use 
the  same  IOT's  for  the  magnetic  tape  device. 

The  processing  performed  on  the  airborne  computer  is 
divided  into  two  parts:  the  first  is  the  collection  of  data 
while  the  flight  is  taking  place,  and  the  monitoring  of  that 
data.  This  first  part  is  called  the  real  time  analysis. 

The  second  use  of  the  airborne  computer  is  referred  to  as 
post  flight  processing  and  is  really  a data  verification 
printout . 
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6.1  DEC  PDP8/E  OPERATION 

This  section  describes  the  procedures  used  to  bootstrap  the 
computer  and  start  up  the  OS/8  operating  system.  OS/8  is  used 
to  call  out  programs  for  execution  such  as,  RTS8 , QWIKY,  etc. 

The  procedures  for  mounting  and  dismounting  magnetic  tapes  are 
also  described. 

Bootstraping  the  Computer 

press  down  HALT  and  pull  back  up 

Dial  in  one  DECtape  unit  as  'O'  (zero),  the  other  as  ' 1' 

REMOTE/OFF/LOCAL  - OFF 

WRITE  ENABLE/WRITE  LOCK  - WRITE  LOCK 

Mount  current  system  tape  on  unit  zero 

Thread  tape  onto  take  up  reel 

REMOTE/OFF/LOCAL  - LOCAL 

press  until  tape  is  secure  on  take  up  reel 
REMOTE/OFF/ LOCAL  - REMOTE  (REMOTE  SELECT  Light  should  come  on) 
set  SR  (switch  register)  to  7470 

press  LOAD  ADDRESS,  EXTD  ADDR  LOAD,  CLEAR,  CONTINUE 
the  hardware  ROM  will  now  cause  OS/8  to  be  loaded  into  core 
and  started 

the  computer  should  print  a dot  (.)  on  the  Tektronix  to  signi- 
fy that  the  OS/8  is  running  and  is  ready  to  accept  commands 

Mounting  and  Dismounting  Magnetic  Tapes 

The  first  step  in  mounting  a tape  is  deciding  whether  or  not 
a write  enable  ring  should  be  used.  The  write  enable  ring 
should  be  used  if  the  tape  is  being  mounted  prior  to  running, 
RTS/8  for  a data  collection  mission.  The  ring  should  not  be 
used  when  running  QWIKY,  IBOL  or  any  program  which  retrieves 
data  already  on  the  tape. 
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Next,  the  following  steps  should  be  followed  to  mount  the 
tape 

turn  transport  on 
LOAD/BR  REL  - LOAD 
ON  LINE/OFF  LINE  - OFF  LINE 
START/STOP  - STOP 
F WD/ RE W/ REV  - FWD 
START/STOP  - START 

the  tape  should  now  advance  to  the  load  point  and  the  LOAD 
POINT  light  should  come  on 
UNIT  SELECT  - 0 

ON  LINE/OFF  LINE  - ON  LINE  (SEL  and  RDY  lights  should  come  on) 
the  tape  is  now  ready  for  use 

To  dismount  tapes  use  the  following  procedure 


ON  LINE/OFF  LINE  - OFF  LINE 
START/STOP  - STOP 
FWD/REW/ REV  - REW 
START/ STOP  - START 

tape  should  rewind  and  position  itself  at  the  load  point 
(LOAD  POINT  light  will  come  on) 

LOAD/BR  REL  - BR  REL 

the  tape  should  now  be  spun  off  by  hand 


6.2  REAL  TIME  SYSTEM  RTS/8 


The  real  time  data  collection  system  consists  of  a monitor 
(RTS/8)  capable  of  managing  a multi-task  environment  and  several 
tasks  which  run  under  control  of  this  monitor.  The  operator  may 
selectively  fun  and  halt  these  tasks  as  desired.  The  primary 
task  which  runs  under  RTS/8  is  the  TAPE  program  which  records  all 
pertinent  data  onto  magnetic  tape.  The  usefulness  of  the  RTS/8 
monitor  lies  in  the  fact  that  it  allows  programs  to  run  real  time 
in  a multi-task  environment,  that  is,  tasks  may  respond  to  exter- 
nal stimuli  and  run  at  the  same  time. 


The  RTS/8  monitor  is  the  heart  of  the  data  collection  system 
Programs  TAPE  and  LWCD  are  run  under  control  of  this  monitor 
To  run  RTS/8  the  following  procedure  should  be  followed 

TEKTRONIX  R4010 
POWER  - ON 
LOCAL/LINE  - LINE 

PAGE  - press  after  screen  lights  up 
HARD  COPIER  4610 
POWER  - ON 

Select  a magnetic  tape  with  the  write  enable  ring  in  to  be 
used  to  record  the  flight  data 
Mount  the  tape  on  the  TU10  (refer  to  section  1.2  for  a 
description  of  this  procedure) 

Bootstrap  the  computer  (refer  to  section  1.1  for  a des- 
cription of  this  procedure)  and  respond  to  the  OS/8 
dot  by  typing  .R  RTS 8 

The  system  should  print  an  asterisk  (*)  on  the  Tektronix 
to  signify  that  RTS/8  is  running  and  ready  to  accept 
commands 
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RTS/8  will  accept  the  following  commands  at  the  command 
level  (signified  by  the  asterisk  printout) 

ON, name  activate  user  program  "name" 

OFF, name  de-activate  user  program  "name" 

To  get  to  command  level  strike  a character  v/hile  RTS/8  is 
running 

To  terminate  data  collection  operations  push  the  HALT  switch 
on  the  computer  and  dismount  the  magnetic  tape  from  the  TU10 
(refer  to  section  1.2  for  a description  of  this  procedure) 

6.2.1  Program  TAPE 

The  TAPE  program  is  the  routine  which  actually  writes  the 
flight  data  out  to  magnetic  tape.  To  begin  the  data  col- 
lection respond  to  the  RTS/8  asterisk  by  typing 

*0N , TAPE 

The  tape  will  then  begin  recording  data.  A record  will  be 
written  every  three  (3)  seconds. 

6.2.2  Program  LWCD 

The  LWCD  program  provides  the  capability  of  monitoring  liquid 
water  content  in  real  time.  Liquid  water  content  is  calcu- 
lated, averaged,  and  displayed  for  a time  interval  defined  by 
the  operator.  Five  different  particle  types  are  supported 
and  the  cloud  or  precip  probes  may  be  used  exclusively  if 
desired . 

To  run  LWCD  respond  to  the  RTS/8  system  asterisk  by  typing 
*ON , LWCD 


: 
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LWCD  will  then  print 
VELOCITY  (KNOTS) 

respond  by  typing  the  velocity  in  knots  (an  integer)  and 
return.  Note  that  all  inputs  to  LWCD  are  terminated  by 
return . 

LWCD  then  prints 

SAMPLE  PERIOD  (SECS) 

respond  by  typing  the  sampling  period  to  use.  This  must 
be  a multiple  of  3.  (3,6,9,12,  etc.)  During  a run  of  LWCD 

the  average  liquid  water  content  during  each  sampling  in- 
terval is  displayed  on  the  Tektronix  screen.  After  the 
sample  period  has  been  entered  LWCD  will  print 

PARTICLE  TYPE 


respond  by  typing  the  particle  type  number  (an  integer  from 
1 to  5) . The  particle  types  currently  supported  by  LWCD 
are: 


1 Rain 

2 Small  snow 

3 Large  snow 

4 Wet  snow 

5 Bullet-Rosettes 

LWCD  then  prints 
PROBE 

respond  by  typing  C to  use  only  the  cloud  probe 

P to  use  only  the  precip  probe 
or  just  return  to  use  both 


-oct; 


in 


LWCD  is  now  running  and  will  display  liquid  water  content 
for  every  sample  period. 

After  every  sample  period  the  following  line  of  output  is 
displayed  on  the  screen. 

HH:MM: SS  TYPE  X LWC  X . XX 


where  HH:MM:SS  is  the  flight  time  clock 

TYPE  X is  the  particle  type  number  being  used 
in  liquid  water  content  calculation 
(input  by  operator) 

LWC  X . XX  is  the  average  liquid  water  content  in 
grams  per  cubic  meter  over  the  sampling 
period  input  by  the  operator 

To  terminate  a LWCD  run  respond  to  the  system  asterisk  by 
typing 

*OFF , LWCD 

(To  get  the  system  asterisk  strike  any  key  while  LWCD  is 
running . ) 
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6.3  REAL  TIME  SYSTEM  RTX/8 

In  addition  to  the  support  and  maintenance  DPSI  has 
provided  on  the  initial  system  (RTS/8),  we  have  designed 
and  developed  a completely  new  real  time  operating  system 
RTX/8.  As  the  year  progressed,  improvements  that  would 
enhance  both  the  quality  and  quantity  of  the  data  collected 
became  more  apparent.  After  considerable  discussion  and 
planning  the  consensus  of  opinion  was  to  create  a new  sys- 
tem rather  than  completely  revise  the  original.  Externally, 
to  the  user,  both  systems  appear  similar;  however  the  level 
of  sophistication  utilized  in  the  new  software  provides  for 
far  more  versatility  as  future  requirements  change,  as  well 
as  an  improved  data  collection  system. 


The  heart  of  the  new  RTX/8  system  utilizes  the  concept 
of  event  processing.  An  event  is  a two-state  condition 
(either  true  or  false)  that  is  used  to  keep  track  of  a par- 
ticular occurence  while  the  system  is  running.  (The  reader 
should  keep  in  mind  that  many  of  these  occurences  take  place 
in  a given  instance.)  The  generality  of  this  event  scheme 
is  extremely  useful  in  a real  time  environment  where  many 
independent  functions  occur  simultaneously.  When  these 
functions  reach  a significant  point  in  their  execution,  the 
computer  must  recognize  this  fact  and  initiate  whatever  action 
is  required.  As  an  example  of  a significant  point  in  the 
processing  stream,  consider  the  task  of  filling  a data  buf- 
fer. As  the  data  is  being  collected  by  the  system  and 
stored  in  the  buffer,  a significant  point  occurs  when  the 
data  buffer  becomes  full.  At  this  time  the  "full  data  buf- 
fer" event  flag  is  set  to  true.  Setting  this  flag  true  sig- 
nals another  program  to  record  the  buffer  on  magnetic  tape. 


After  the  data  buffer  is  recorded  the  event  flag  is  cleared 
(i.e.  set  to  false)  and  the  process  is  repeated.  The  fol- 
lowing list  shows  some  of  the  on  going  functions  that  uti- 
lize this  event  flag  technique. 

1.  PMS  synch  pulse  detected 

2.  Tektronix  screen  erased 

3.  PMS  buffer  processing  completed 

4 . Operator  typed  1 C on  keyboard 

5.  Liquid  water  content  calculation  completed 

6.  Magnetic  tape  error  detected 


The  versatility  of  this  technique  is  apparent  from  this 
list  of  functions.  Tasks  to  the  system  may  be  added  or  de- 
leted easily  without  affecting  the  overall  processing  flow. 
This  event-flag  scheme  has  another  powerful  advantage  over 
the  earlier  system.  In  certain  applications  the  status  of 
one  event  is  important  to  many  different  programs  and  spe- 
cialty routines. 


Consider  the  operation  of  the  Tektronix  screen.  This 
device  is  used  by  many  different  programs  (plus  the  opera- 
tor) during  data  collection.  When  the  screen  becomes  full 
it  must  be  copied  and  erased;  all  programs  using  the  screen 
must  be  made  aware  of  that  fact.  The  plotting  program,  for 
example,  must  redraw  the  axes  and  begin  plotting  at  the 
origin.  Only  one  event  flag  is  required  to  inform  all  the 
other  user  programs  that  the  screen  has  been  erased. 


The  event-flag  technique  allows  many  functions  to  be 
performed,  that  would  have  been  impossible  otherwise.  Some 
of  these  capabilities  are  listed  below; 
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I 

Simultaneous  usage  of  both  the  Tektronix  screen 
and  GE  printer  as  independent  input/output  devices 

Detecting  and  reporting  all  errors  due  to  magtape 
problems  or  data  transmission 

Reporting  all  inoperable  devices 

Monitor  integrity  - i.e.  user  programs  cannot 
cause  system  crashes 

Complete  communication  between  all  user  programs 
if  desired 

Full  exploitation  of  this  concept  offers  numerous  advantages 
that  were  unavailable  to  the  earlier  system.  £t  this  time 
RTX/8  is  fully  operational;  additional  features  will  be  in- 
stalled in  the  near  future. 


/ 


I 

I 


j 
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6.4  POST  MISSION  DATA  REDUCTION  SYSTEM 

The  post  mission  data  reduction  system  consists  of  a series 
of  programs  which  read  ID  PMS  tapes  or  PDP8  generated  tapes  and 
perform  various  lookup  and  display  options.  Specific  time  in- 
tervals may  be  retrieved.  These  programs  would  be  used  on  site 
for  preliminary  data  analysis  prior  to  processing  on  the  AFGL 
CDC  6600. 


6.4.1  Proram  QWIKY 


QWIKY  performs  a quick  look  type  dump  for  all  ID  PMS  and 
PDP8  data  tapes.  PMS  tapes  may  be  used  from  the  C130A  or 
the  C130E.  The  dump  printed  by  the  program  includes  time 
printouts,  raw  and  calibrated  VCO's,  total  probe  counts, 
and  certain  values  derived  from  this  data. 


The  program  may  be  run  in  two  different  modes.  The  first 
mode  dumps  the  entire  tape  printing  averages  over  a speci- 
fied interval.  The  second  mode  allows  flight  time  or  elapsed 
time  to  be  used  to  locate  a record.  This  feature  will  be 
useful  in  locating  the  start  of  a sampling  run  or  any  time 
interval  of  interest.  The  operating  instructions  for  QWIKY 
are  shown  on  the  following  four  pages. 


Mount  the  magtape  to  be  dumped  on  unit  0 with  the  write 
enable  ring  removed  to  prevent  any  possible  corruption  of 
data 


Turn  on  the  GE  Terminet  printer,  press  the  'ON  LINE'  button 
and  be  sure  that  the  paper  is  free  to  feed. 


Respond  to  the  OS/8  dot  (.)  by  typing  "R  KNMON"  (return). 

/ 

(for  detailed  explanation  of  how  to  get  the  OS/8  dot  refer 
to  section  6.1,  Bootstrapping  the  Computer. 


When  QWIKY  is  started  it  will  print 


QUICK-LOOK  PROGRAM 

FLIGHT  TAPE  (F)  OR  KENNEDY  TAPE  (K) ? 


Respond  by  typing  F if  the  magtape  is  a PDP8  generated  flight 
tape  of  K if  the  tape  is  from  the  Kennedy  tape  recorder  in 
the  PMS  system.  Terminate  input  with  a return.  Note:  all 
inputs  to  QWIKY  are  terminated  with  return.  If  K is  input 
the  program  will  then  print 


Respond  with  A if  the  tape  is  from  the  C130A  or  E if  the 
tape  is  from  the  C130E.  The  program  will  then  print 


FLIGHT  NO 


DATE? 


Enter  the  flight  number  and  date  or  a string  of  up  to  30 
characters  which  identifies  the  tape.  The  program  then  asks 


TIMES? 


Respond  with  A to  dump  the  whole  tape  starting  at  whatever 
record  the  tape  head  is  currently  at.  The  program  will  not 
rewind  the  tape  to  the  load  point  before  dumping.  This  is 
useful  if  the  program  fails  to  locate  a specified  time  in- 
terval. The  operator  may  position  the  tape  to  approximately 
where  the  desired  time  whould  be  using  the  off  line  controls 
on  the  TU10.  The  A option  may  then  be  used  to  determine 
what  was  recorded  on  the  tape.  There  is  one  small  problem 
associated  with  this  procedure.  When  the  off-line  controls 
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are  used  to  forward  space  or  backspace  the  tape  the  TU10 
does  not  stop  on  an  even  record  boundary.  Therefore  the 
first  few  seconds  of  data  produced  by  the  program  should 
be  ignored  if  tape  has  been  moved  with  the  off-line  con- 
trols. The  program  will  detect  this  error  condition  and 
print 


!!! INCORRECT  MAGTAPE  RECORD  LENGTH!!! 

This  printout  is  normal  after  the  tape  has  been  moved  with 
the  off-line  controls,  however  it  should  only  print  once . 
The  first  read  done  by  the  TU10  will  reposition  the  tape 
head  to  the  beginning  of  a record  and  subsequent  records 
should  be  readable  without  generating  this  error.  If  this 
printout  repeats  continuously  the  tape  is  not  in  the  cor- 
rect format  and  any  results  printed  by  the  program  will  be 
incorrect . 

The  S option  allows  the  operator  to  specify  a time  interval 
on  the  tape  to  be  dumped.  If  this  option  is  selected  the 
program  then  prints 

CLOCK  (C)  OR  ELAPSED  (E)  TIME? 


Type  C to  use  the  flight  time  clock  to  look  up  records  or 
E to  use  the  PMS  buffer  elapsed  time.  The  program  will 
attempt  to  find  the  specified  time  interval  no  matter  where 
it  is  on  the  tape  and  will  minimize  tape  motion  as  much  as 
possible.  If  the  tape  has  gone  beyond  the  time  interval 
desired  the  program  will  backspace  as  necessary.  It  is  not 
necessary  to  rewind  the  tape  to  the  load  point  before  each 


After  the  operator  has  specified  which  clock  to  use  the 
program  will  print 


START  TIME 


Respond  by  typing  the  starting  time  for  the  dump  in  the  form 
SSSS  to  specify  a four  digit  elapsed  time  in  seconds  or 
HH : MM: SS  to  specify  a flight  time.  The  program  will  then 
print 


STOP  TIME 


The  program  now  asks 


AVERAGING  INTERVAL  (SECS) 


This  question  is  printed  for  both  the  specified  time  interval 
and  dump  all  times  mode.  Respond  with  the  number  of  seconds 
to  be  represented  in  each  printout. 


The  program  then  begins  calculating  averages  over  the  speci- 
fied time  interval  and  printing  a report  for  each  interval. 
Contained  in  each  printout  are  the  start  and  stop  flight  and 
elapsed  times  for  the  interval,  average  VCO  values,  average 
values  of  parameters  derived  from  the  VCO's,  and  total  probe 
counts . 


The  program  may  be  restarted  from  the  beginning  by  setting 
the  switches  to  0200  and  pressing  HALT,  ADDR  LOAD,  EXT  ADDR 
LOAD,  CLEAR,  CONTINUE. 


Alternatively  the  program  can  be  restarted  from  the  "ALL  (A) 
OR  SPECIFIED  (S)  TIMES?"  question  by  setting  the  switches  to 
0000  and  pressing  HALT,  ADDR  LOAD,  EXT  ADDR  LOAD,  CLEAR, 
CONTINUE. 
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Error  Recovery  Procedures 


The  following  error  messages  may  print  out  during  execution 
of  QWIKY . 


printout 

possible  cause (s) 

recommended  action ( s ) 

TIME  INTERVAL  OUT  OF 
RANGE 

time  specified  not 
on  tape 

use  "A"  option  to  see 
what  is  on  tape 

**END  OF  TAPE** 

time  specified  not 
on  tape 

rewind  tape  and  try 
another  time  interval 

1 ! ! MAGTAPE  PARITY 
ERROR! ! ! 

bad  tape,  bad  TU10 

try  another  tape; if 
the  printout  is  in- 
frequent it  may  be 
ignored /however  if 
data  looks  bad  this 
is  the  cause 

! ! ! INCORRECT  MAGTAPE 
RECORD  LENGTH ! ! ! 

bad  tape  (unlikely) 
question  asking  if 
tape  was  Kennedy 
tape  (PMS)  or 
Flight  tape  (PDP8 ) 
was  answered  in- 
correctly tape  head 
was  not  on  even  re- 
cord boundary  due 
to  operator  moving 
it  using  off-line 
controls 

try  another  tape; re- 
start program  at  0200 
and  answer  question 
correctly 

no  action  necessary 
simply  ignore  results 
of  first  averaging 
interval 

- 

Some  errors  will  not  be  detected  by  the  program.  Following 
is  a list  of  the  conditions  and  recommended  actions. 


disposition 

possible  causes (s) 

recommended  action  (s) 

TUlO  "rocks"  and  will 
not  locate  time  speci- 
fies 

bad  times  on  tape 

use  "A"  option  to 
dump  tape 

program  will  not  load 

bad  system  tape 
unknown 

try  backup  system  tape 
notify  DPSI 

program  "hangs  up" 

unknown 

record  AC,  PC,  MQ  and 
notify  DPSI 

program  halts 

unknown 

record  AC,  PC,  MQ  and 
notify  DPSI 

6.4.2  Program  IBOL 


IBOL  provides  the  means  of  retrieving  flight  data  from  flight 
data  tape  created  by  the  RTS/8  real  time  data  collection  pro- 
gram. Data  from  the  Knollenberg  probes  is  specified  by  VCO 
number  and  may  be  retrieved  on  the  basis  of  time,  referencing 
either  the  aircraft  clock  or  the  elapsed  time  counter  of  the 
Knollenberg  device.  Data  may  be  printed  on  the  line  printer 
or  displayed  on  the  CRT  (DECwriter  for  Lab  version)  as  speci- 
fied by  the  operator. 

To  run  IBOL  mount  the  PDP8  generated  flight  tape  with  the 
write  enable  ring  removed  to  prevent  any  possible  corruption 
of  data 

Turn  on  the  GE  Terminet  printer,  press  the  'ON  LINE'  button 
and  be  sure  that  the  paper  is  free  to  feed. 

Respond  to  the  OS/8  dot  (.)  by  typing  "R  IBOL"  (return). 


Once  initiated  by  the  R IBOL  directive  the  IBOL  program  re- 
quests input  from  the  operator  regarding  the  flight  identi- 
fi ration,  reference  clock,  start  and  stop  times  for  the  data 
desired  and  finally  a list  of  the  Knollenberg  words  to  be 
retrieved  and  displayed. 

Once  the  required  parameters  have  been  supplied,  the  program 
IBOL  will  commence  a search  of  the  records  of  the  RTS/8 
flight  tape  until  the  data  of  the  specified  time  period  is 
located.  The  desired  KNW  values  will  then  be  retrieved  and 
displayed  on  either  the  CRT  or  line  printer  depending  upon 
which  has  been  specified  as  the  output  device  as  follows: 

Set  -.witch  register  bit  zero: 

SR  bit  0 = 0 output  on  CRT  (DECv:riter  for  lab  version) 


SR  bit  0=1  output  on  Line  Printer 

When  all  the  data  for  a given  time  interval  has  been  dis- 
played, the  program  will  initiate  itself  to  allow  a new  pass 
The  data  tape  is  not  rewound  at  this  time,  (see  example  of 
program  operation) 

operator  inputs  (all  underlined  characters  are  program  gene- 
rated) 

A)  FLIGHT  IDENTIFICATION 


NOTE:  MAKE  OUTPUT  DEVICE  SELECTION  NOW 

operator  may  enter  up  to  500  characters  to  provide  an 
identification  text  which  will  appear  on  the  output  (if 
the  line  printer  is  selected  as  the  output  peripheral, 
or  on  the  CRT  if  it  is  the  specified  output  device.) 

Any  character  may  be  entered  with  the  exception  of  the 
back  slash  which  is  used  to  terminate  entry. 

NOTE:  Errors  during  entry  of  the  flight  id  can  be  cor- 
rected by  use  of  the  'rubout'  key.  Each  time  the  rub- 
out  key  is  depressed,  the  last  character  is  erased  from 
the  identification  string.  To  restart  the  entire  id 
entry  operation,  simply  depress  'control'  and  ' C'  simul- 
taneously . 

B)  ENTER  AC  OR  ET  FOR  CLOCK  SELECTION 

operator  responds  by  typing  AC  (space)  if  the  aircraft 
clock  is  to  be  used  as  the  time  reference  for  data  inter 
val  specification  ET  (space)  if  the  elapsed  time  from 
Knollenberg  is  to  be  used  as  the  time  reference. 

C ) BEGIN  . . . 

enter  start  time  of  the  interval  over  which  data  is  to 
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be  retrieved.  Two  formats  are  used  depending  upon  the 
clock  selection  made  in  step  (B)  . 

1)  if  AC  was  selected,  the  program  expects  a time  given 
in  the  form  of  hours,  minutes  and  seconds: 

HH : MM: SS  (space) 

2)  if  ET  was  selected,  the  program  wants  an  elapsed  time 
in  seconds  between  0 and  9999: 

1234  (space) 


D)  END 


enter  the  finish  time  of  the  interval  over  which  data  is 
to  be  retrieved.  The  form  of  time  entry  must  be  the  same 
as  step  (C) . 

E)  ENTER  KNW  SELECTION  FOR  DISPLAY 
TERMINATE  OPERATION  WITH  A SPACE 


operator  may  enter  the  KNW  numbers  to  be  retrieved  and 
displayed.  The  range  of  KNW  numbers  is  1 to  63.  Selec- 
tion of  KNW  numbers  outside  that  range  will  result  in  in- 
correct values. 

KNW  numbers  may  be  specified  as  individual  numbers  sepa- 
rated by  commas  or  as  a sequence  of  contiguous  KNW  ' s by 
means  of  a dash  to  separate  the  first  from  the  last  KNW. 
In  each  case,  the  KNW  number  or  sequence  must  be  followed 
by  a comma  with  the  exception  of  the  last  number  which  is 
followed  by  a slash  to  terminate  the  operation. 

example:  To  display  KNW  numbers  1,6  thru  12,9,22  and 
45  thru  63,  the  following  type  in  is  made: 
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NOTE: 

Errors  made  during  the  entry  of  KNW  numbers  may  be  cor- 
rected by  depressing  the  'rubout'  key.  This  has  the 
effect  of  clearing  any  partially  input  number  and  allow- 
ing a new  entry.  Once  the  number  has  been  completed  how- 
ever, (terminated  with  a comma,  space  etc.)  the  rubout 
will  have  no  effect.  Note  also  that  operation  of  the 
rubout  key  will  not  result  in  carriage  motion  or  charac- 
ter echo. 

Once  the  required  information  has  been  provided  to  the 
program,  the  desired  data  is  located  on  the  RTS/8  flight 
tape  and  transmitted  to  the  device  selected  as  the  out- 
put peripheral  (CRT  or  line  printer) . 

Data  is  recovered  until  the  end  of  the  interval  speci- 
fied in  step  D is  reached.  At  such  time,  the  message 
'END  DATA'  is  typed  and  the  program  automatically  re- 
initiates itself  for  another  run.  If  for  some  reason, 
the  run  must  be  interrupted  and  new  parameters  used, 
striking  any  key  will  cause  the  program  to  re-initialize 
itself  after  the  current  display  operation. 


On  completion  of  any  given  run,  the  flight  data  tape 
is  positioned  at  the  point  of  last  data  acquisition.  If 
it  is  necessary  to  recover  other  KNW  values  from  an  earlier 
time  period  on  the  same  flight  tape,  the  tape  must  be 
manually  rewound  prior  to  commencing  the  run. 
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6.5  TESTING  and  CALIBRATION 


j 


t 


i 

j 

> 

' 
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Shortly  before  a scheduled  data  collection  mission,  it 
is  desirous,  whenever  possible,  to  perform  certain  system 
tests.  This  insures  that  the  data  collected  will  be  of  the 
best  quality  prossible.  Many  minor  system  malfunctions  can 
be  corrected  at  this  time.  Also,  periodically  throughout 
the  season,  the  VCO  calibrations  are  rechecked.  The  program 
KNMON  was  written  to  assist  the  flight  crew  when  performing 
these  tests. 
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6.5.1  Program  KNMON 

This  program  was  developed  for  testing  and  checking  out 
of  the  PMS  interface  and  associated  M1703  card  in  the  compu- 
ter. Any  desired  Knollenberg  channel  may  be  selected  for 
monitoring  by  the  operator.  The  channel  selected  is  printed 
every  second  since  the  PMS  system  sends  its  blocks  of  infor- 
mation at  this  rate.  Since  the  channel  is  printed  every  time 
the  PMS  system  sends  it  to  the  computer,  printout  from  the 
program  should  be  useful  for  calibration  work  and  trouble^ 
shooting  VCO ' s in  addition  to  checkout  of  the  basic  PMS  buf- 
fer functions  and  the  sync  word. 

To  operate  the  program 

1.  Respond  to  the  OS/8  dot  by  typing  "R  KNMON"  (return) . 

2.  After  KNMON  is  started  it  will  print 

KN  CHANNEL 

Respond  by  typing  the  Knollenberg  channel  to  monitor 
(from  1 to  64)  followed  by  return.  The  program  will 
then  print  the  value  of  the  selected  channel  every 
second . 

3.  To  select  another  channel  set  SR  to  0200  and  press  'HALT', 
'LOAD  ADDRESS',  ' EXTD  ADDR  LOAD',  'CLEAR',  'CONTINUE'  then 
continue  from  step  2. 


______ 
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7 . ADDITIONAL  TASKS 


During  the  year  certain  tasks  became  apparent  that 
would  ease  the  everyday  processing,  and  in  some  cases  the 
analysis,  that  occurs  at  LYC.  Usually  these  take  the 
form  of  short  real-time  programs  to  be  used  on  the  CDC  6600 
terminal;  however,  in  some  cases  the  DEC  PDP8/I  and  E are 
utilized.  In  all  cases  the  goal  is  to  eliminate  some  of 
the  tedious  but  necessary  work  required  of  the  LYC  scientists. 


I 
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7.1  Program  LOOP 

Program  LOOP  is  the  last  analysis  performed  in  the  post 
flight  processing.  As  a result  of  the  runs  of  KN0LL1D  and 
RAPP,  the  relationships  between  radar  reflectivity  and  liq- 
uid water  content  are  calculated  by  fitting  functions.  A 
missile  will  be  fired  into  this  same  cloud  formation  after 
the  aircraft  has  landed,  and  radar  observations  will  be 
taken  as  the  missile  passes  through  the  atmosphere.  The 
radar  response  will  be  values  of  reflectivity  for  which 
liquid  water  content  will  be  desired.  (These  z values  are 
taken  from  the  LYC  "storm  top"  plot.) 

The  input  to  program  LOOP  is  handl  d by  a special  pro- 
gram called  BOBIN . Both  LOOP  and  BOBIN  are  designed  to  be 
executed  at  the  6600  terminal,  with  output  printing  at  the 
computer  center  • 


BOBIN  will  first  request  the  21  points  which  describe 
number  density  versus  size.  These  21  numbers  are  output 
from  KNPLTlD . The  diameters  of  these  21  points  are  fixed 
at  .075,  .15,  .25,  and  from  .4  to  3.8  spaced  at  .2.  The 
program  will  also  request  the  three  sets  of  A,  B coeffi- 
cients from  KNPLtID  and  RAPP.  These  three  A's  and  B's  are 
described  below: 

* 
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B3/3 

D o = A3(Za/Ma)  from  KNPLT1D 

where  the  subscript  'a'  indicates  an  aircraft  derived 
parameter  and  the  subscript  ' r'  indicates  a radar  derived 
parameter 


Program  BOBIN  allows  the  user  to  print  and  modify  any 
of  these  27  values  (21  points  from  KNPLTlD,  2 sets  of  A,  B 
from  RAPP  and  A3,  B3 . When  satisfied,  control  is  tranferred 
to  program  LOOP  for  execution.  This  separation  of  function 
was  done  in  order  to  allow  LOOP  a restart  function  without 
having  to  reinput  the  27  numbers.  It  is  especially  useful 
for  rerun,  and  in  the  event  System.  L at  the  computer  center 
crashes.  Program  BOBIN  creates  a permanent  disk  file  of  these 
numbers;  it  also  purges  the  previous  file  when  a new  file  is 
made  and  updates  the  cycle  number.  This  data  file  is  called 
X4340. 

Program  LOOP  will  determine  the  mass  (M)  and  median 
volume  diameter  (D0)  for  measured  radar  reflectivity  using 
empirical  distributions  of  dimensionless  diameter  (D/D0) 
vs.  dimensionless  number  density  (NDq^/M) . A short  scien- 
tific paper  has  been  written  describing  the  techniques  for 
this  determination.  This  paper  in  its  entirety  is  included 
here . 
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From  the  running  of  analysis  programs  on  flight  data  (see 
KN0LL1D  and  RAPP  in  "DEVELOPMENT  AND  APPLICATION  OF  DATA  PROCES- 
SING TECHNIQUES  AND  ANALYTIC  PROCEDURES  TO  CLOUD  PHYSICS  DATA", 
by  Belsky , Kaplan  and  Rodenhiser,  AFCRL-TR-7 5-04 27 , 28  July  1975) 
three  fitting  functions  are  produced.  The  first  yields  a fit  of 
22  points  of  D/Dq  (equivalent  melted  diameter  over  median  volume 
diameter)  by  ND04/M  (number  times  the  fourth  power  of  median  vol- 
ume diameter  divided  by  the  liquid  water  content) . The  second 
produces  a least  square  fit  of  the  parameter  (mk)  such  that 

(mk)  = Ai  ZBl 

where  (mk)  is  defined  as  the  liquid  water  content  divided  by  the 
square  root  of  the  radar  reflectivity,  Z is  the  radar  reflectivity 
and  Aj  and  Bx  are  the  fitting  parameters.  The  third  produces  a 
least  square  fit  to  the  median  volume  diameter  such  that 

B2/3 

D0  = A2  (M/Z) 

where  M is,  again,  the  liquid  water  content  and  the  A2  and  B2  are 
fitting  parameters. 


Our  task  is  for  a given  (D/Do)_3_,  and  a given  radar  reflec- 

IUaX 

tivity,  Z , to  find  the  D0  and  M which  satisfy  the  equations 


ZR=  / D 


6 _M_ 
DO4 


FdD 


■ifl 


M = z I D6  FdD 

6 n ** 

0 D0 


Zdr.t  = /(D/D°)maX  D6  — — — FdD 


MOBJ  6 ^ 


(D/D  o ) 


D3  _EL_  FdD 


MOBJ  (lak)  ’/ZOBJ 


F is  the  output  function  ND04/M 

M is  the  mass,  or  liquid  water  content 

D0  is  the  median  volume  diameter 

Z is  a measured  radar  reflectivity 
.K 

D is  an  integration  variable 

(D/Dn)  is  the  instrument  maximum  (will  be  referred 
' u max 

to  as  IM) 


The  reported  output  for  F,  or  ND04/M  is  actually  their  logs 
base  10.  It  will  thus  be  more  cc  .ivenient  to  write  lO-*7  for  F, 
where  y's  are  the  numbers  actually  being  printed.  In  addition, 
the  third  fitting  function: 


Dn  = A,  Md2 


Bp/3 


will  be  used  solely  to  obtain  an  original  estimate  to  D0 


The  analysis  for  solution  follows. 
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given  the  following  equations 


M 


M 

2 


M 

2 


7T 

6 


/ D3  10y  dD 

0 D04 


= l fD°  D3  SL.  I0y  dD 


M , v 


6 1 
0 


D, 


= \ I D3  — 10y  dD 


= / D6  10Y 


dD 


MOBJ  = (mk) /ZOBJ 


IM 

mobj  - ? / °3  rr  loy  aD 

0 Do 


IM 


OBJ 


/ D6  -A-  10Y  dD 
0 D04 


investigating  the  first  three  equations,  it  is  helpful  to  replace 
10y  with 


-kD/D ( 


Nn  e 


(see  "Snow  Size  Spectra  and  Radar  Reflectivity"  by  Sekhon  & 
Srivastava  Journal  of  American  Sciences,  Vol  27,  No  2 March 
1970  pp  299-307) 


W 
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then 


—iuo  I D3e-kD/D»  dD 

Do4  6 J 


„ / D3e-kD/D»  dD 

Do'  6 0 


« , JL  ” N„  / Dle-W>°  dD 

D»“  D 

Do 


dividing  2A  by  3A  yields 


/DVe-kD/D» 


/ D3e‘kD/D0  dD 


this  can  be  integrated  directly,  yielding 


6D04  D04  3D  0 4 6D04  6D  0 4 

- ( + + + ) e 

k4 k k^ k_^ k4 

D04  3Dq  4 6D04  6D0  4 

( + + + ) e 

k k2  k3  k4 


this  can  be  simplified  to 


and  solved  for  k.  In  this  case  k = 3.67206075 


The  more  important  result  is  that  maintaining  equation  4A  is  in- 
dependent of  D o . That  is,  provided  D0  is  chosen  finite  and  greater 
than  zero,  this  relationship  must  be  maintained.  This  suggests 
that  our  function  be  changed  from 


where  a is  chosen  such  that  the  fol 


lowing  equation  holds 


This  equation  is  also  independent  of  the  choice  of  D0*  This 
remarkable  fact  can  be  seen  from  adding  equations  2A  and  3A, 
yielding 


again  choosing  D0  = 1. 


From  equations  1 to  7 we  have 


(Iili_)2  < 
6D04 


/ 


D3  10ay  dD  f 


(mk)2  — — / d6  10ay  dD 

D04  0 


and  solving  for  M 


(mk) 2 

M = Sal 

©’  (; 


IM 

/ D6  I0ay  dD 


IM 

/ D3  10aY  dD 


I 


or 


M 


IM 

(mk)2  / D6  10ay  dD 

0 


D n 4 36 


IM 

/ D3  10“y 


dD 


combining  this  result  with  equation  4B,  we  have  to  solve  for  Do 
such  that 


Z 


R 


(Jl_)  e / d6  10ay  dD 
Do4  0 


(mk)  2 


IM 

/ D6 


10aY  dD 


where 
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The  system  yields 

a,  B which  are  independent  of  Z , and  depend 

R 

only  upon  the  given  distribution  Y(D/D0) 

Do  the  solution  of  the  preceding  equation 
M an  auxiliary  value  in  the  preceding  equation 
as  defined  by  7B 
Mqbj  as  defined  by  5 


138 


1.  by  Newton's  method,  find  a such  that 


/ D3  10ay  dD 

0 


= / D3  10ay  dD 


by  forming 


F(a)  = / D3  10aY  dD  - / D3  10ay  dD 

0 1 


then  letting  a = a0f  a starting  value  (generally  1) 


form 


F(ak)  , F(1.01ak) 


then 


F'  = 100[F(1.01ak)  - F(ak)]/ak 


and 


ak+l  = ak  “ F(ak)/F' 


the  process  is  continued  until  ak+-^  - ak  -*■  0 


2.  Find  3 


ir  / D3  10  Y dD 
0 


3.  By  Newton's  method,  find  D0  such  that 


then  letting  (D0)0  = A2  M 


then 


process  is  continued  until  (D0) 


Find  M 
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The  operating  instructions  for  these  two  programs  are 
found  on  the  following  pages. 

To  enter,  alter,  print  data  file  (BOBIN) 

COMMAND-ATTACH , DATA,  BOBINX4  34  0,  ID=BELSKY* 

'if  LFN  ALREADY  IN  USE,  or 

PF  ATTACHED  go  to  next  step 

COMMAND-REQUEST, NFILE , *PF* 

if  get  duplicate  message:  COMMAND-UNLOAD, NFILE, FILE* 

COMMAND- REQUEST, NFILE, *PF* 

COMMAND-DATA 

if  program  doesn't  run:  COMMAND- UNLOAD, DATA*! 

go  to  top  line  J 


To  execute  calculating  program  (LOOP) 

COMMAND-ATTACH , LGO , LOOPX  4 34  0 , ID^BELSKY* 
if  LFN  ALREADY  IN  USE,  or 

PF  ATTACHED  go  to  next  step 


COMMAND-LGO 

if  CANNOT  ATTACH  DATA  FILE:  COMMAND-UNLOAD, NFILE, FILE 

COMMAND-LGO 

To  end  run,  type  Zr=0.  (ALWAYS  TYPE  DEC.  PT. 

IN  RESPONSE  TO  ZR, INSTR.MAX,MK) 

To  get  rid  of  output:  COMMAND-DISPOSE, TAPE6* 

After  running  LOOP:  COMMAND-UNLOAD, NFILE, FILE* 

To  Rerun  BOBIN:  COMMAND-REQUEST , NFILE , *PF* 

COMMAND-DATA* 

To  Rerun  LOOP:  COMMAND-LGO* 


A 
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H 


to  erase  a character  type  CTRL/H;  one  character  is  erased 
for  each  time  the  keys  are  depressed, 
to  erase  an  entire  line,  type  CTRL/X;  then  type  as  if 
nothing  had  been  previously  typed. 


LOGIN:  LOGIN,  LOREN,  I.YB44 9 , SUP-t 

LOGIN,  MORT,  LIAl  7 8 , SUP  4- 
LOGIN, FRED, L20461, SUP+ 


trouble  X4546  (operator)  X4005  (D.  Caveliero) 


IF,  in  BOBIN,  the  message 

CANNOT  CATALOG  DATA  FILE.  ERROR  NO.  015 
all  is  not  lost!  simply  type 
COMMAND -CATALOG, NFILE,X4 3 4 0, ID=BELSKY , CY=1 + 
COMMAND-PURGE, FILE+ 

COMMAND-UNLOAD, FILE! 


142 


1 


7 . 2 Program  F.AIN 

RAIN  is  the  program  for  processing  the  raindrop  spectro- 
meter data.  This  data  is  recorded  in  analog  form  from  the 
raindrop  spectrometer  directly.  The  analog  tape  is  run  on 
the  LYC  PDP8I  computer  via  the  analog  to  digital  converter, 
and  the  digitized  results  are  written  on  the  magnetic  tape 
and  printed  on  the  high  speed  printer. 

The  output  from  the  PDP81  is  run  by  program  RAIN  on  the 
CDC  6600  in  order  to  produce  the  desired  results.  The  format 
for  the  records  on  the  magnetic  tape  will  be  found  in  appenc ix 
38.  These  records  are  preceded  by  a "header"  record  of  52 
alphabetic  characters,  and  they  are  followed  by  an  end-of-file. 
The  26  channels  contain  counts  of  raindrops  for  1 second,  by 
diameter  class. 

The  diameter  class  depends  upon  the  number  of  bit:  (k) 

used  in  the  analog  to  digital  conversion.  This  value  of  k 
will  have  to  be  furnished  to  program  RAIN,  and  is  given  by 
the  A to  D operator.  The  operator  is  working  with  voltage 
values  on  the  analog  tape.  He  will  determine  k in  such  a 
way  that  the  voltage  values  will  be  divided  into  the  26 
classes  without  obtaining  values  beyond  this  range.  At  the 
same  time,  if  k is  too  large,  all  values  will  be  in  the  be- 
ginning of  the  26  classes  without  sufficient  separation. 

That  is,  too  large  a k yields  too  small  a resolution. 

The  program,  itself,  will  calibrate  these  results,  calcu- 
late liquid  water  content  (N),  reflectivity  (Z)  and  rain 
rate  (R) . It  will  also  plot 
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M vs . Z 
R vs . Z 
M vs . R 
time  vs.  M 
time  vs.  R 
time  vs . Z 

Finally  the  last  plot  will  show,  for  each  P second  period, 
number  density  (N)  vs.  diameter,  with  a Marshall-Palmer  line 
fit  through  the  spectrum. 

I 

To  determine  the  diameters  for  each  of  the  26  channels: 

d = /(22 . 5)  m + 4.3  (microns) 

m 2^ 

where 

n = channel  number 

k = number  of  bits  in  A to  D conversion 


and  we  define  do  = /4.3 

The  midpoint  diameter  value  (D^)  for  any  one  of  the  26  classes 
is  given  by 

d . + d 

D = -2 — 51  (microns) 

m 2 

It  is  interesting  to  note,  as  an  aside,  that  the  maxi- 
mum diameter  is  a function  of  k,  the  number  of  bits  in  the 
A to  D conversion.  The  minimum  diameter  is  always  2.0736 
microns  (dQ),  and  the  maxima  for  some  selected  k values  are: 
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maximum  diameter  (microns) 


7.0713 


5.2107 

3.9656 

3.1643 


The  formulas  used  in  the  calculation  are: 


Number  density,  N 


count 


N = 
m 


(counts/m3 ' 


where  v is  sampling  volume 


Liquid  water  content,  or  mass,  M 


fl  = \ N o D 3 (gms/m3' 
m 6 m m 


where  p = water  density 


Radar  reflectivity,  Z 


Z = N D 6 
m mm 


(mm6/m3 ) 


Rain  rate,  R 


R = 
m 


(3600)  (14.2)  / m 


— (mm/hr ) 


1000 
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and  these  are  summed  over  the  26  classes  to  produce  results 
for  the  1 second  interval: 


1 = I M 

L i 
m 

- I Z 

m 

- I 


There  are  two  peculiarities  of  the  magnetic  tape  that 
cause  difficulty.  The  first  is  that  channel  1 gives  spurious 
results  and  the  data  for  the  first  channel  must  be  ignored. 
Secondly,  the  time  on  the  tape,  (hours,  minutes  and  seconds) 
is  often  incorrect.  This  second  fault  requires  that  time 
will  have  to  be  controlled  externally. 

Since  the  time  cannot  be  relied  upon,  the  user  will 
supply  the  time  of  the  first  data  record,  and  each  subsequent 
record  will  have  a time  of  one  second  later,  until  an  end-of- 
file  is  seen.  Thus,  there  will  have  to  be  a beginning- time 
given  for  each  file.  These  times  can  be  determined  from  the 
high  speed  printout  of  the  PDP8I  during  the  writing  of  the 
magnetic  tape. 


7 . 3 Program  HELP 


HELP  is  a program  designed  to  run  in  a conversational 
mode  on  the  6600  time  sharing  system  INTERCOM.  The  program 
will  produce  a punched  deck  or  disk  file  containing  all  the 
input  required  to  run  the  programs  KNOLL1D  or  KNPLTlD . The 
output  consists  of:  a job  card,  control  cards,  and  data  cards. 
It  is  particularly  useful  when  many  decks  have  to  be  generated. 
A disk  file  is  produced  that  allows  the  job  to  be  submitted 
through  the  BATCH  remote  entry  system.  Operation  of  this  pro- 
gram is  quite  simple;  the  sample  run  on  the  following  pages 
contain  all  the  information  necessary  for  execution. 


' 
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COMMAND-  ATTACH  HELP,  HELPS  I NX2876,  I D=GLASS 
COMMAND-  CONNECT, INPUT 

COMMAND-  CONNECT, OUTPUT 
COMMAND-  HELP 

ENTER  DESIRED  PROGRAM  NAME ...  KNOLL  1 D 

JOB  NAME (5  CHAR) 7GLAA8 

(B) INART  OR  (S)OURCE?  B 

CYCLE, 1 3 OR  BLANK  FOR  LATEST.. 209 

PROGRAM  NAME. .KNOLL! D 

INPUT  TAPE  NO.  ? LYC136 

PLOT  TAPE  NO.?  LYC178 

RAPP  TAPE  NO.?  LYC179 

CHANGING  DCOEF?  N 

CHANGING  PCOEX?  N 

CHANGING  SCOEF ? N 

( A ) OR  <E)  FLIGHT?  A 

DO  YO'J  WISH  TO  SEE  THE  VCOEF.?  N 

CHANGING  VCOEF?  N 

START  TIME-  ENTER  AS  HHMMSS • • 2 1 0900 

STOP  TIME-  ENTER  AS  HHMMSS  ..  22 1 300 

FLIGHT  ID  ENTER  AS  XYR-NN . . A7 6- 1 A 

FLIGHT  DATE-  ENTER  AS  DD  MON  YR..25  APR  76 

PM S OMTIME  AS  HHMMSS . . 2 1 0 1 57 

TYPING  BY  (H)ANUAL  INPUT  OR  <T)EMP..M 

(A)/C  CLOCK  OR  ( P >MS  BUFFER.. A 

(P  ) T'.r’.L  l <■  OR  ( F ) I i\'AL  DATA  . . F 

V;riAT  V E INTERVAL  (IN  SECONDS  < 12))  • .(Vi 


TEMP  BY  (V)CO^  C A ) T MO SPHERE.  OR  (R)ADIOSONDE 

CALCULATE  UC 0?  YES 

SKIP  KOw  M AMY  FILES  (12).. 30 

PUNCH  HECK?  NO 

(K)WAJ  OR  ( S ) r ANDARD?  S 

DATA  MODIFICATI ON ?N0 

<T)YPF  OR  (L)DIT  OR  (N)EITHER. . T 

START  TIME-  ENTER  AS  HHMMSS . . 2 I 1 000 

STOP  TIME-  ENTER  AS  HHMMSS •• 2 1 1 530 

FOR  CLOUD  CODE. .07 

FOR  PRECIP  CODE.. 07 

(T)YPE  OR  (E)DIT  OR  ( N ) El THER • • T 

START  TIME-  ENTER  AS  HHMMSS . .21 1 600 

STOP  TIME-  ENTER  AS  HHMMSS ..  2 1 1 800 

FOR  CLOUD  CODE.. 07 

FOR  PRECIP  CODE.. 07 

(T)YPE  OR  (E)DIT  OR  ( N ) El THER . . T 

START  TIME-  ENTER  AS  HHMMSS •• 21 2200 

STOP  TIME-  ENTER  AS  HHMMSS . .2 1 2900 

FOR  CLOUD  CODE. .07 

FOR  PRECIP  CODE.. 07 

CT)YPF  OR  (E)DIT  OR  (N)EITHER. . T 

START  TIME-  ENTER.  AS  HHMMSS  ••  2 1 5800 

STOP  TIME-  ENTER  AS  HHMMSS . .220200 

FOR  CLOUD  CODE..  3 -5 

FOR  Pn'.CIP  CODE.. 05 

( T ) Y P OR  (F)DIT  "R  ( N ) FI  THEN . .N 


COMMAND-  ATTACH  HELP*  HELPS I NX 2 8 76* I D=  GLASS 
COMMAND-  CONNECT/ INPUT 
COMMAND-  CONNECT/ OUTPUT 
COMMAND-  HELP 

ENTER  DESIRED  PROGRAM  NAME. . .KNFLT1 D 

</».) A AM  '••  ( S C R LAP 9 

CYCLE ( 1 3)/  OP.  ->L ALA  FOR  LATEST... 

INPUT  TAPE. . .LYCI7P 

(B) INARY  OR  ( S ) OURCE.  . .B 

(L>  EAP.JET  OR  (O130...C 

( P ) EN  OR  (ORT..C 

(A)/C  OR  (P)MS  CLOCK... A 

CA)LL  DATA  OR  JUST  ( S ) UMMARY . . . S 

FLIGHT  ID  ENTER  AS  XYR-NN. . . A76- 1 A 

FLIGHT  DATE  AS  DD  MON  YR..25  APR  76 

JLAR  SET  OR  (D)EFINE  WHICH  ONES  TO  PLOT . . K 

THIS  LOOP  WILL  ASK  FOR  EACH  PASS  INFORMATION  WHEN  LAST  PASS  IS  GIVEN  TYPE 
-1  FOR  PASS  NO. 

THEN  THE  HEIGHT  PROFILES  WILL  DE  GENERATED  FROM  THE  PREVIOUS  INFOP_MATION 
PASS  NO.< 12) 01 

START  PASS  TIME  AS  HHMMSS  • . ft  1 1 600 
STOP  PASS  TIME  AS  HHMMSS .. 2 1 1 800 
HEIGHT (KILO)  TIMES  TEN  ( T • D -TD  ) IN  ( I 3 ) • • 089 
PASS  J0-CI2) 0 P. 

START  PASS  TIME  AS  HHMMSS .. 2 1 2200 
STOP  PASS  TIME  AS  HHMMSS .. 2 1 2900 
HEIGHT (KILO)  TIMES  TEC  ( T • D"TD)  I?  (13). .076 
»AS;  .(12). ..21 

ST  A;  i ( f I HP  •i;K 5S  • » R 1 5830 

stop  ms  T i • - as  .22. 

i,  • ( f LO)  TIMES  TEN  (T.D-TO)  IN  (1.3). .OS? 


7 . 4 Program  PLOT 


Program  PLOT  was  written  for  problems  arising  at  LYC, 
where  the  data  was  not  recorded  on  magnetic  tape.  Also, 
LYC  had  requirements  for  a plotting  program  which  were 
dissimilar  to  the  requirements  for  the  on  board  computer. 


Program  PLOT  will 

(a)  plot  an  x-y  table  on  the  Tektronix  plotter;  the 
table  is  inoutted  at  the  Decwriter.  If  desired, 
the  x and/or  y values  can  be  modified  by  logging 
them  in  order  to  produce  a linear-linear,  a linear- 
log,  a log-linear  or  a log-log  plot. 

(b)  generate  the  x values  automatically  after  a speci- 
fication of  the  first  x and  x-step  is  given 

(c)  allow  the  user  to  select  the  low  and  high  values 
of  x and  y to  be  used  on  the  plot 

(d)  plot  each  point  with  +;  the  user  can  choose  whether 
or  not  to  connect  the  plotted  points  with  a line 

(e)  label  the  plot  with  an  alphabetic  description 

(f)  plot  a least  square  best  fit  curve  to  the  data 
(first  or  second  degree)  , and  print  out  the 
fitting  function 

(g)  generate  a table  of  deviations  of  the  least  square 
fit  to  the  original  data,  and  calculate  the  RMS 
error 

(h)  allow  the  user  to  modify  the  data,  limits  and 
descriptions  and  replot,  with  a new  least  square 
fit,  without  having  to  retype  the  entire  x-y 
table . 
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In  order  to  execute  program  PLOT,  and  perform  the  many 
optional  capabilities,  the  following  step-by-step  procedure 
should  be  executed: 

1.  Mount  Dectape  136  on  unit  0,  write  enabled  and  remote 

2.  Switches  set  to  7‘.70 

3.  Press  Addr-load,  Ext . -addr-load,  Clear,  Cont 

4.  Computer  responds  with  a dot 

5 . Type  R PLOT 

6.  Within  22  seconds  the  computer  responds  with 

PLOTTING  PROGRAM 

7.  IS  X ON  LOG  SCALE? 

If  the  x-data  is  to  be  logged  base  10,  answer  Y otherwise  N 

8.  IS  Y ON  LOG  SCALE? 

If  the  y-data  is  to  be  logged  base  10,  answer  Y otherwise  N 

9.  IS  DELTA  X CONSTANT? 

If  the  x-data  (before  taking  logs)  is  equally  spaced, 
answer  Y,  otherwise  N 

10.  If  delta  x is  constant,  the  computer  will  request  X-START 
and  DELTA-X,  as  X-START? 

Type  the  starting  value  and  return, 

DELTA-X? 

Type  the  difference  between  successive  x-values  and  return 

11.  The  computer  responds  with 

ENTER  X, Y TABLE 

or  ENTER  Y TABLE  (for  constant  delta-x) 

■ 


i 
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12.  For  each  value  input  the  computer  will  first  print  the 
index  number  ( 1 , 2 , 3 , . . . etc ) If  only  y-values  are  being 
inputted,  the  user  will  type  the  y-value  followed  by  a 
return;  if  both  x and  y values  are  being  inputted,  the 
user  will  type  the  x-value,  followed  by  a comma,  followed 
by  the  y-value  and  a return. 

13.  When  all  values  have  been  inputted,  type 

1E35  return 

14.  The  computer  will  respond  with 

LIMITS  OF  X ARE  xxxxx,  xxxxx 
which  indicate  the  low  and  high  values  of  the  x-table. 
Then  it  prints 

TYPE  LIMITS  TO  USE 

The  computer  expects  two  numbers,  separated  by  a comma , 
the  first  number  is  the  value  of  x at  the  left  end  of 
the  plot;  the  second  is  the  value  of  x at  the  right  end 
of  the  plot.  There  will  automatically  be  10  divisions 
in  this  range.  If  the  first  limit  typed  is  unacceptable 
(larger  than  the  limit  found)  or  if  the  second  limit 
typed  (smaller  than  the  limit  found)  is  unacceptable, 
step  14  is  repeated. 

15.  Step  14  is  done  for  Y 

16.  The  computer  will  print 

CONNECTING  LINE? 

If  a line  is  desired  between  points  in  the  order  they 
were  typed  in,  answer  Y.  An  answer  of  N will  eliminate 
tile  connecting  line. 

17.  The  computer  will  respond 

TURN  ON  PLOTTER,  TYPE  DESCRIPTION 
Type  a description  to  be  printed  on  the  bottom  of  the 


plot.  This  description  is  limited  to  67  characters.  A 
rubout  will  respond  with  a carriage  return,  and  the  user 
should  type  the  entire  description  over  again. 

Make  sure  plotter  is  on,  and  "on-line"  before  pressing 
return . 

13.  The  plot  will  be  made.  At  its  conclusion,  the  computer 
will  ask 

LEAST  SQUARE  FIT? 

If  a curve  fitting  the  data  is  desired  answer  Y ; other- 
wise answer  N and  proceed  to  step  22. 

19.  The  computer  responds 

DEGREE (1  OR  2)? 

Degree  1 is  a straight  line  fit;  degree  2 is  a quadratic 
fit.  Respond  either  1 or  2 . Any  other  response  will 
force  step  19  again. 

20.  The  fitting  curve  will  be  plotted.  The  least  square  poly 
nomial  will  be  printed 

(a)  at  the  top  of  the  plot,  and 

(b)  on  the  Decwriter 

21.  The  computer  then  asks 

DEVIATIONS? 

If  the  table  of  x,  y,  calculated-y , deviation-y  and  RMS 
is  desired,  answer  Y;  otherwise  answer  N 

22.  The  computer  will  ask 

REPLOT? 

If  the  data  is  to  be  modified,  eliminated,  or  extended, 
or  if  another  plot  is  to  be  made  changing  either  x or  y 
limits,  or  the  description,  or  the  fitting  function, 
respond  with  Y;  otherwise  respond  N and  proceed  at  step  6 
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23.  The  computer  will  ask 
NUMBER? 

If  there  are  no  changes  go  to  step  23(c). 

(a)  The  computer  expects  the  index  number  of  the  data 
point  being  changed.  It  will  accept  an  index  num- 
ber one  greater  than  the  length  of  the  table  (omit 
the  number  corresponding  to  1E35)  with  the  assump- 
tion that  the  table  is  being  extended.  The  table 
may  be  extended  repeatedly  by  entering  an  index 
number  equal  to  one  more  than  the  previous  maxi- 
mum number.  If  the  number  typed  is  invalid,  step 
23  is  repeated.  After  the  number  is  accepted,  the 
computer  will  type  X,Y  = 

Respond  with  the  value  of  x,  a comma,  the  value  of 
y and  return.  Step  23  will  be  repeated.  The  data 
being  typed  will  be  logged  according  to  responses 
made  in  steps  7,8 

CAUTION:  DO  NOT  ENTER  1E35  FOR  THE  LAST  NUMBER. 

The  program  knows  how  many  points  there  are  in  the 
table . 

(b)  If  a value  of  x,y  is  to  be  eliminated,  after  its 
index  number  has  been  accepted,  and  x,Y  = has  printed, 
type 

1E35,1E35  return 

(c)  When  all  changes  have  been  made,  answer  NUMBER?  with 
0 return 

The  procedure  restarts  at  step  14 


Ifc-J* 


ERRORS 

1.  Alphabetic  errors:  If  an  error  is  discovered  before  re- 

turn is  pressed,  and  the  information  was  alphabetic,  as 
in  steps  7, 8, 9, 16, 17, 18 , 19 , 2] , 22,  type  the  rubout  key; 
the  computer  will  respond  with  a carriage:  return;  retype 
the  entire  response,  beginning  with  the  first  character. 

2.  Numeric  errors:  If  an  erro;  is  discovered  before  a comma 

or  return  is  pressed,  and  the  information  was  a number, 
as  in  steps  10,12,14  and  23,  type  the  rubout  key;  the 
computer  will  not  respond  (automatic  feature  of  the  27 
bit  floating  point  package);  then  type  the  entire  number 
over  again,  beginning  with  the  first  character.  CAUTION: 
If  tv/o  numbers  were  being  inputted,  separated  by  a comma, 
the  rubout  will  erase  the  CURRENT  NUMBER  ONLY. 

3.  Incorrect  responses  to  Y or  N:  The  computer  will  speci- 

fically look  for  the  Y.  Any  character  other  than  Y will 
be  treated  as  N except  a carriage  return,  which  should 
be  used  only  after  a character  has  been  typed. 

RESTART  ADDRESS:  0200 


Operating  Instructions 


The  deck  setups  for  the  various  programs  follow 
sections  8.1  through  8.7: 


H 

00 

KNOLL ID 

CN 

00 

KNPLTlD 

ro 

00 

KNlUTIL 

CO 

HIAC1D 

8.5 

SPANDAR 

00 

CT\ 

RAPP 

8.7 

RAIN 

i 


8.1  KN0LL1D 
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PROGRAM  KNOLL ID 


DECK  SETUP  (OBJECT  FILE) 


CONTROL  CARDS 


JOBNR, CM110000 , T900 , TPN . * Prob.  No.  Name 

ATTACK , LGO , KNOLL IDE I NX287  6 , ID-GLASS , MR= 1 . 

ATTACH , A1 0 , Al 0 FORMAT . 

REQUEST  TAPE1 , S , H I , VSN=PKSXXX . 

REQUEST  TAPE 2 , RING, VSN~TAPENO . 

REQUEST  TAPE?, RING, VSN=TAPENO. 

FILE  ( TAPE1 , RT=U, BT=K, MRL=1 0 2 4 , MBL=1 0 24 , RB=1 , BFS=1 05) 
LDSET, FILES=TAPE1 , PRESET^ZERO . 


MAP, OFF. 


1,0 AD , Al  0 . 

LGO. 

REWIND , TAPE 3 , TAPE 8 . 

COPY, TAPE 3 . 

COPY , TAPE8 . 

EXIT. 

REWIND, TAPE 3 , TAPE8 . 

COPY, TAPE3 . 

COPY , TAPE8 . 

7/8/9 

DATA  CARDS 

7/8/9 

G/7/8/9 

= Nun, o<  r of  REQUEST  TAPI  cards  (i.c.  1,3,  or  3) 
for  PLottapo  option  must  have  REQUEST  TAPF.2  card 
for  Rapptape  option  must  have  REQUEST  TAPE?  card 


l 


•| 
: j 


I 
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DATA  CARD  SEQUENCE  SUMMARY 

A.  REQUIRED  DATA  CARDS 
card  # 1 s 
1-4 

5 

6 

7 

8 

3.  OPTIONAL  DATA  CARDS 


(9-18) 

JVCO  cards 

(OPTIONAI.-only 

if 

card 

8 

max 

cc  34-35  is  greater  than 

zero) 

(19-38) 

RADIOSONDE 

(OPTIONAL- only 

if 

card 

8 

max 

cards 

cc  30  is  two) 

(39-43) 

EDIT  cards 

(OPTIONAL) 

max 

(44-58) 

TYPE  cards 

(OPTIONAL) 

max 

TYPE 

and  EDIT  cards  may 

be  intermixed 

coefficient/exponent  cards 
start /stop  processing  time 
tape  card 

flight  information  card 
option  card 


h. 


. . . I y ~ 
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COEFFICIENT/EXPONENT  CARDS 

The  four  cards  $DCOEF,  $PCOEX,  5SC0EF  and  $VCOEF  are 
used  to  change:  the  adjusted  class  coefficients  (with  their 
breakpoints) ; equivalent  melted  diameter  coefficients  and 
exponents  (with  their  breakpoints);  sounding  coefficients; 
and  VCO  coefficients. 

The  breakpoints  allow  the  coefficients  or  exponents 
for  a given  particle  type  to  be  changed  as  the  channel  num- 
ber or  crystal  size  gets  larger.  BRKPTN  controls  the  coef- 
ficients of  the  adjusted  class  equation;  the  units  are  in 
N (channel  number).  BRKPTC  controls  the  coefficients  and 
exponents  of  the  equivalent  melted  diameter  equation;  the 
units  are  in  mm  (crystal  size) . 

The  particle  type  code  is  shown  as  follows: 

J Type 

RAIN 
WET  SNOW 
LARGE  SNOW 
SMALL  SNOW 
BULLET-ROSETTES 
COLUMNS 
NEEDLES 

RJMC  EXPERIMENT  51 
RMC  EXPERIMENT  #2 
UNUSED  AT  THIS  TIME 


1 

3 

5 

7 

9 

11 

13 

15 

17 

19 


1.  $DCOEF  card 

This  card  allows  the  coefficients  of  the  adjusted  class 
equation 

N'  = mN  + b where  N = measured  class 

to  be  changed.  The  breakpoints  for  this  equation  are  also 
changed  using  this  card. 

Control  variables:  K BRKPTN ( J ) 

where:  I = 1 or  2 and  J = 1 to  9 

1=1  for  the  slope  (n) 

1=2  for  the  intercept  (b) 

J = particle  type  code  (N  < breakpoint) 

J+l = particle  type  code  (N  > breakpoint) 

For  example,  to  set  the  slope  for  type  rain  equal  to  1.0 
for  classes  1 through  3 and  equal  to  1.5  for  classes  4 
through  15  use: 

@ $DC0EF  MB  ( 1 , 1 ) = 1 . 0 , MB ( 1 , 2 ) = 1 . 5 , BRKPTN ( 1 ) = 3 , BRKPT ( 2 ) = 3 $ END 

NOTE:  This  breakpoint  is  probe  independent  i.e.  the  break- 
point occurs  twice,  once  for  the  cloud  probe  and 
once  for  the  precip  probe. 
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2.  $PCOEX  card 

This  card  allows  the  coefficients  and  exponents  of  the 
equivalent  melted  diameter  equation 

0X 

D = co (d  ) where  d = crystal  size  in  mm 

to  be  changed.  The  breakpoints  for  this  equation  are  also 
changed  using  this  card. 

Control  variables:  COEX ( I, J) , BRKPTC { J) 

where  I = 1 or  2 and  J = 1 to  9 

I =1  for  the  coefficient  (co) 

I - 2 for  the  exponent  (ex) 

J = particle  type  code  (d  < breakpoint) 

J+l  = particle  type  code  (d  > breakpoint) 

For  example,  to  set  the  exponent  for  type  Small  Snow 
equal  to  1.5  for  crystal  sizes  less  than  or  equal  to  .4  mm 
and  equal  to  2.0  for  crystals  sizes  greater  than  .4  mm 
use: 

@ $PC0EX  COEX ( 2 , 7 ) =1 . 5 , COEX ( 2 , 8 ) =2 . 0 , BRKPTC ( 7 ) -BRKPTC ( 8 ) = 0 . 4 $ END 

NOTE:  This  breakpoint  is  probe  dependent  i.e.  the  break- 
point occurs  only  once  regardless  of  the  probe 
being  used. 
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3.  $SCOEF  card 

This  card  allows  the  coefficients  of  the  pressure-height 
sounding  equation  to  be  changed. 


Control  variable:  S(I) 


where  I = 1 to  5 


first  order  coefficient 
second  order  coefficient 
third  order  coefficient 
fourth  order  coefficient 
fifth  order  coefficient 


For  example,  to  change  the  second  order  coefficient, 
the  control  variable  is  S(2)  = new  value. 

$VCOEF  card 

This  card  allows  the  VCO  calibration  coefficients  to  be 
input,  xhe  default  coefficients  are  0,1  and  0.  Thus 
each  run  must  have  the  desired  coefficients  input. 

Control  variable:  C(I,J) 

where  I = 1 to  3 and  J = 1,13  for  A Model 

J = 1,7  for  E Model 


A MOD 


intercept  i pressure 

slope  2 delta-p 

third  order  coef.  3 heading 

4 temperature 

5 event  and  cloud 

6 LWC-JW 

7 rain  spectrometer 

8 Tacan  bearing 

9 Tacan  distance 

10  acceleration 

11  dewpoint  (1011) 

12  icing  rate  meter 

13  pitch 

For  example,  to  input  the  third  order 
he  control  variable  is  C(3,6)  = value. 


E MOD 

delta-p 
temperature 
pressure  (msb) 
pressure  (lsb) 
dewpoint  (1011) 
LWC-JW 
heading 


coefficient  of  JW, 


w 


k 
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5.  Start  & stop  processing  card 

These  variables  control  the  data  to  be  processed;  any 
time  prior  to  the  start  time  or  after  the  stop  time  is  ig- 
nored . 


cc 

1-6 

START= 

cc 

7-16 

SSI? 

cc 

17-20 

blank 

cc 

21-25 

ST0P= 

cc 

i 

26-35 

SS@ 

6.  Tape  card 

This 

i 

card  contains  the  input  and  output 

If  desired 

it  may  be 

left  blank. 

cc 

1-6 

INPUT  TAPE  NUMBER 

cc 

7-10 

blank 

cc 

11-16 

PLOT  TAPE  NUMBER 

cc 

17-20 

blank 

cc 

21-26 

RAPP  TAPE  NUMBER 

1 


I 
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7.  Flig 


cc 


cc 


cc 


cc 


it  information  card 


1-10  Flight  Id  FLT0XYR-NN 

where:  X must  be  aircraft  model 

i . e . A or  E 

YR  is  the  last  2 digits  of  the 
year 

NN  is  the  flight  number 


11-20  Flight  Date  @DD@MON@YR 

where:  DD  is  the  day  of  the  month 

MON  is  the  first  three  letters 
of  the  month 

YR  is  the  last  2 digits  of  the 
year 


21-24  blank 


25-30 


PMS  on  Time  HHMMSS 

where:  HHMMSS  is  the  time  corresponding 


to  PMS  count  of  0000 


r 

I 

K 


8.  Option  card 

This  card  specifies  the  options  to  be  used.  It  also 
causes  certain  optional  cards  to  be  required. 


J cc 

5 

ICRYS 

= 

1 

TYPING  BY  TEMPERATURE 

= 

2 

UNUSED 

1 

= 

3 

TYPING  BY  MANUAL  INPUT 

cc 

10 

ICLK 

1 

TIME  FROM  AIRCRAFT  CLOCK 

= 

2 

TIME  FROM  PMS  BUFFER 

i cc 

15 

IDAT 

= 

0 

PRELIMINARY  DATA 

> 

0 

FINAL  DATA 

1 

I cc 

20 

IPLT 

2= 

0 

PLOT  TAPE  NOT  PRODUCED 

1 

SS 

1 

PLOT  TAPE  PRODUCED 

cc 

22-25 

INTA 

2S 

0 

FULL  DATA 

> 

0 

AVERAGING  INTERVAL  IN  SECONDS 

cc 

30 

ITMP 

2= 

0 

TEMPERATURE  DETERMINATION  BY  VCO 

= 

1 

DETERMINATION  BY  STANDARD  ATMOSPHERE 

2 

TEMPERATURE  DETERMINATION  BY 
RADIOSONDE  PROFILE 

1 

cc 

35 

JVCO 

= 

0 

CALCULATE  VCO  PARAMETERS 

> 

0 

NUMBER  OF  JVCO  CARDS 

cc 

40 

NSKP 

0 

PROCESS  FIRST  FILE 

> 

0 

NUMBER  OF  FILES  TO  SKIP  BEFORE 
PROCESSING 

cc 

45 

IDEK 

= 

0 

DECK  NOT  PUNCHED 

' 

= 

1 

DECK  PUNCHED 

cc 

50 

IRAD 

= 

0 

RAPP  TAPE  NOT  PRODUCED 

= 

1 

RAPP  TAPE  PRODUCED 

cc 

55 

KWAJ 

- 

0 

STANDARD  PROBE  CONFIGURATION 

• 

= 

1 

KWAJALEIN  PROBE  CONFIGURATION 

Option  card  (cont'd) 


NOTES 


ICP.YS 

For  ICRYS=3  at  least  1 TYPE  card  is  required 
ICLK 

For  ICLK=1  PMS  on  time  may  be  set  to  zero 
IDAT 

For  IDAT>0  IDAT  is  printed  as  FINAL  £ on  page  2 of 
output. 


ITMP 

For  ITMP=1  Standard  atmosphere  equation  is: 

T = 76.88288 (p)0-150284 

For  1 TMP=2  Profile  uses  a maximum  of  20  levels 


JVCO 

For  JVCO> 0 allows  VCO  parameters  to  be  manually  input 
for  specified  time.  This  option  overrides  the  ITMP 
selection 
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9-18  (max)  JVCO  cards  (optional) 

JVCO  cards  are  required  only  if  cc  34-35  on  the  option 
card  is  greater  than  zero.  These  cards  override  the  VCO 
calculations;  they  supply  the  basic  meteorological  parameters 
necessary  for  processing.  The  input  values  are  constant  for 
the  duration  of  the  time  specified.  The  card  format  follows 


cc  •• 

FORMAT 

DESCRIPTION 

1-10 

A10 

@HH:MM:SS@  interval 

start  time 

11-20 

sAlO 

@HH: MM: SS@  interval 

stop  time 

21-30 

. F10.0 

PRESSURE  (mb) 

free  field 

31-40 

F10. 0 

AIRSPEED  (knots)  ' i 
TEMPERATURE  (°C)  ‘ 

1 j 

with 

41-50 

F10. 0 

decimal 

51-60 

F10.0 

HEIGHT  (feet)  j L 

point 
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19-38  (max)  RADIOSONDE  cards  (optional) 


These  cards  are  required  only  if  the  option  card  cc  30 
is  2.  They  override  the  VCO  temperature  calculation  by 
allowing  the  temperature  to  be  determined  by  3 inear  .inter- 
polation using  the  RADIOSONDE  temperature-pressure  profile. 
The  first  card  contains  the  number  of  cards  to  follow  (20 
maximum) . The  remaining  cards  contain  the  RADIOSONDE  pro- 
file. These  cards  must  be  in  descending  pressure  sequence. 


card  1 


cards  2 . . . n 


cc  1-2 

cc  1-10 
cc  11-20 


number  of  levels 


PRESSURE  (mb) 
TEMPERATURE  (°C) 


39-43  (max) 


EDIT  cards  (optional) 


These  cards  allow  selected  particle  channel  counts  to 
be  edited.  The  current  editing  constraints  are  as  follows 
To  edit  channel  1 or  15  - there  must  be  two  valid 

channels  on  one  side 

To  edit  channels  2-14  - there  must  be  one  valid 

channel  on  both  sides 
To  edit  channels  1 and  2-  there  must  be  a valid 

channel  3 


cc 

FORMAT 

DESCRIPTION 

1-4 

A4 

EDIT 

5-14 

A10 

START  TIME 

(@HR:MN: SC?) 

15-24 

A10 

STOP  TIME 

(@HR:MN:  SC>?) 

25-26 

12 

PROBE  NO . 

f 1 = scatter 

27-23 

12 

PROBE  NO. 

| 

> < 2 = cloud 

29-30 

12 

PROBE  NO.  J 

| 

v.  3 = precip 

31-33 

13 

CHANNEL  NO 

34-36 

13 

CHANNEL  NO 

37-39 

13 

CHANNEL  NO 

40-42 

13 

CHANNEL  NO 

43-45 

13 

CHANNEL  NO 

46-48 

13 

CHANNEL  NO 

49-51 

13 

CHANNEL  NO 

• 

52-54 

1 3 

CHANNEL  NO 

. 

L 


44-58  (max) 


TYPE  cards  (optional) 


The  type  cards  are  required  only  if  the  option  card 
cc  5 is  3.  These  cards  specify  the  particle  type  to  be 
used  for  the  selected  interval.  The  cloud  probe  type 
and  precip  probe  type  may  be  different.  The  scatter  probe 
is  always  type  rain. 

TYPE  CARDS 


cc 

FORMAT 

DESCRIPTION 

1-4 

A4 

TYPE 

5-14 

A10 

START  TIME  ( @HR: MN : SC@ ) 

15-24 

A10 

STOP  TIME  ( @HR: MN : SC@ ) 

25-26 

12 

TYPE  KEY  (Cloud  Probe) 

27-28 

12 

TYPE  KEY  (Precip  Probe) 

TYPE  KEY 

1 

RAIN 

3 

WET  SNOW 

5 

LARGE  SNOW 

7 

SMALL  SNOW 

9 

BULLET-ROSETTES 

11 

COLUMNS  (4:1) 

13 

NEEDLES  (7.5:1) 

15 

RMC  1 -EXPERIMENTAL 

17 

RMC  2 -EXPERIMENTAL 

19 

UNUSED  AT  THIS  TIME 

NOTE:  rp’ne  TYPE/EDIT  cards  are  interchangeable  and  can  be  in 
^.ny  order.  There  is  a maximum  of  15  TYPE  cards  and  5 
EDIT  cards,  any  more  will  be  ignored. 
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KNOLLlD  V2.01 

1.  For  upper  altitude  pressures  (P  < 440mb) 

the  following  pressure  correction  is  applied  ... 

PCOR  = 1 • 25*Pi-110 

where 

PCOR  = corrected  Pressure 
= calculated  Pressure 

2.  For  a/c  time  between  12:07:08  & 12:48:53  use 

AP  = 37 . 06mb  FIXED 

NOTE  - this  version  is  for  flight  A75-15  only (15  APR  75) 

OPERATING  CHANGES 

1.  control  cards: 

use  cycle  = 201 

2.  data  cards; 


none 


n 


^ i 
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KNOLLlD  v2 . 02 

1.  Height  calculated  from  interpolation  of  time-height  profile 

2.  Pressure  calculated  from  ... 

P = P0(l-(|  )HT)1/aR 

where 


Pq  = Pressure  @ ground 
Tq  = Temperature  @ ground 
a = Lapse  rate  @ ground 
R = Gas  constant 
HT  = Height 

OPERATING  CHANGES 


1.  control  cards: 

use  cycle  = 202 

2.  data  cards: 

insert  immediately  after 

1st  card:  c£ 

1-2  (rj) 
11-20 (ff) 
21-30  (ff ) 
31-40  (ff) 
41-50 (ff ) 


option  card 

desc~ - tion 
# of  levels  to  follow 


a 


R 


rj  = right  justified 

ff  = free  field  with  decimal  point 

remaining  cards: 

(50  max)  cc  description 

1-10  time@HH:MM:SS@ 

11-20  (ff)  height  (meters) 

NOTE: 

1.  start  time  must  be  greater  than  or  equal  to  the  time 
of  level  1 

2.  stop  time  must  be  less  than  or  equal  to  the  time  of 
the  last  level 


i 


- 


V 


( 
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KNOLL ID  v2 . 03 

1.  For  Precip  Probe,  type  2 only,  the  breakpoint  occurs  at 

N > 5 instead  of  N > 3.  Other  changes  for  this  probe  are 
done  with  data  cards. 

2.  The  Cloud  Probe  changes  are  also  done  with  data  cards. 

NOTE  - This  version  was  created  from  2.01.  It  should  be 
used  with  flight  A75-19  only  (15  APR  75) 

OPERATING  CHANGES 

1.  control  cards; 
use  cycle  203 

2.  data  cards 
@$DCOEF 

@MB( 1,2) =0.4 7 5, MB (2, 2) =3. 53, 

@MB (1, 3) =1.170, MB (2, 3) =0. 0, 

@ MB  ( 1 , 6 ) = 1 . 1 7 0 , MB  ( 2 , 6 ) = 0 . 0 , 

0BRKPT ( 2 ) =1 . 7 55 , BRKPT ( 3 ) =1 . 755 , 

@ $END 
@ $PCOEX 

0COEX (1, 2)=0. 370,COEX(2, 2) =0. 670, 

0COEX (1, 3)=0. 370, COEX ( 2 , 3)=0. 670, 

(3$  END 
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KN0LL1D  v2 . 04 


r 


I 


1.  This  version  incorporates  the  latest  changes  requested  by 
Dr.  R.  M.  Cunningham  through  24  NOV  75. 

2.  The  particle  type  list  has  been  expanded.  Refer  to  the 
operating  instructions  (page  8)  for  a complete  list  of 
particle  types. 

3.  The  breakpoint  option  has  been  modified  to  include  two 
types  of  breakpoints:  crystal  size  and  channel  number. 
Pages  5 and  5.1  of  the  operating  instructions  reflect 
these  changes. 

NOTE  - This  version  was  created  from  2.03.  It  should  be 
used  v/ith  flight  A75-19  only. 


OPERATING  CHANGES 

1.  control  cards: 
use  cycle  = 204 

2.  data  cards: 


L 


i 


$DCOEF  see  operating  instructions  for  card  1 
TYPE  CARD  see  new  particle  type  key 


KN0LL1D  v2 . 0 5 


1.  This  version  is  the  standard  for  all  flights  before  Decern 
ber  1975,  except  15  APR  75.  It  incorporates  all  the 
changes  of  v2.04. 

OPERATING  CHANGES 

1.  control  cards: 
use  cycle  205 

2.  data  cards: 
none 
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KNOLL ID  v2 . 06 


1.  This  version  is  the  standard  for  all  flights  after 
December  1975. 


2.  It  should  be  used  with  all  PMS  tapes  with  the  new  format. 
The  tape  format  changes  are  summarized  below: 


WORD# 


NEW 


OLD 


1 

17 

33 

49 


elapsed  time 
BO  status  word 
CO  status  word 
DO  status  word 


sync 


probe  configuration 
elapsed  time 
unused 


OPERATING  CHANGES 


1.  control  changes: 
use  cycle  = 206 

2.  data  cards: 
none 


m 


i 
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KN0LL1D  V2.07 

1.  This  version  contains  changes  to  data  modification 
routine  used  in  v2.06.  These  changes,  primarily  for 
processing  the  data  from  flight  A76-04  on  11  February 
1976.,  include  the  following: 

A.  Data  modification  is  not  used  for  particle 
types  RAIN  or  WET  SNOW. 

B.  If  there  are  insufficient  points  from  the  precip 
probe  (i.e.  less  than  two  channels)  the  data 
modification  technique  is  not  used.  However, 
the  cloud-precip  gap  is  still  eliminated  and  M, 

Z,  Dq  and  K are  recalculated.  The  gap  elimina- 
tion technique  has  not  been  changed.  Refer  to 
the  KN0LL1D  Operating  Manual  for  a complete 
description  of  this  technique. 


, 


!* 

J 


OPERATING  CHANGES 

1.  control  cards: 
use  cycle  = 207 

2.  data  cards: 


none 


KNOLLlD  v2 . 08 


1.  This  version  utilizes  the  data  modification  routine 
from  v2.07  and  the  pressure  correction  routine  from 
v2 . 02  . 

2.  The  height  is  determined  by  interpolation  of  the  time 
height  profile. 

3.  The  pressure  is  calculated  from: 

P = P0  (1-  (=£)  HT)  1/aR 
U l0 

where: 

Pq  = pressure  @ ground 
Tq  = temperature  @ ground 
a = lapse  rate  @ ground 
R = gas  constant 
HT  = height 

OPEARATING  CHANGES 

1.  control  cards: 
use  cycle  --  208 

2.  data  cards: 

insert  immediately  after  option  card 

1st  card:  cc 

1-2  (rj) 

11-20 (ff) 

21-30  (ff ) 

31-40  (ff ) 

41-50 (ff ) 


description 
# of  levels  to  follow 


a 

R 


rj  = right  justified 

ff  = free  field  with  decimal  point 


remaining  cards 
(50  max) 


description 


time  @HH:MM:SS@ 
height  (meters) 


NOTE 


start  time  must  be  greater  than  or  equal  to 
the  time  ot  level  1 


stop  time  must  be  less  than  or  equal  to  the 
time  of  the  last  level 
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KNOLL1D  v2 . 09 


1.  This  version  will  produce  a line  printer  plot  of  total 
Dq , M,  and  Z as  a function  of  time.  The  three  variables 
are  plotted  on  the  same  logarithmic  scale  with  a range 
of  + 4.0 

2.  The  output  deck  format  has  been  changed  to  the  following: 


| 

i 


cc 

variable 

1 

blank 

2-7 

time  (HHMMSS) 

8-10 

blank 

11-15 

height  (feet) 

21-30 

LWC  scatter 1 

31-40 

LWC  cloud  V gm/M3 

41-50 

LWC  precip  J 

51-60 

D0  cloud  ^ 

61-70 

D0  precip  v microns 

71-80 

Do  total  J 

3.  The  timing  and  central  memory  requirements  have  been  re- 
duced considerably.  This  should  improve  program  turn 
around  time. 


OPERATING  C RANGES 

1.  control  cards: 

a.  job  card  requires  CM61000  and  a time  limit 
of  150  seconds  per  hour  of  data. 

b.  use  cycle  209 

c.  for  the  line  printer  plot  add  the  following 

cards  after  the  LGO. 

REWIND, TAPE9 . 

COPY , TAPE9 . 

2.  data  cards 

option  card  col  65  = 1 to  produce  the  line  printer 
plot 
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KNPLT1D 


CONTROL  CARDS 

JOBNR,  CHI 6 0 0 0 0 , T50 0 , TP  1 . Prob.  No.  Name 

REQUEST  TAPE2 , VSN=  . 1 

ATTACH, LGO, KNPLT1DBINX2876 , ID=GLASS, CY= ,MR=1 . 

ATTACH, CRT , OFFLINECRT, MR=1 . 

DISPOSE, TAPE39, *FL. 

LIBRARY, CRT, 

LDSET, PRESET-ZERO . 

LGO. 

7/8/9 

DATA  CARDS 

7/8/9 

6/7/8/9 

THE  LEARJET  TAPE  IS  TO  BE  PLOTTED  THIS  CARD  MUST  BE 


REQUEST  TAPEl , VSN= 
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DATA  CARDS 

The  first  data  card  is  the  information  card;  it  appears 
only  once.  The  format  for  this  card  is  shown  below. 

Each  plot  type  requires  a plot  request  card.  These 
request  cards  are  unlimited  and  have  the  same  format.  The 
cards  are  divided  into  16  fields.  Each  plot  type  requires 
certain  fields  to  be  used;  all  the  unused  fields  may  be  left 
blank . 

card  1.  information  card 


VAR 

CC 

FORMAT 

FUNCTION 

PLT 

1-3 

A3 

plot  type:  PEN  or  CRT 

CLK 

5 

11 

clock:  1=A/C  2=PMS 

IOUT 

7 

11 

0 = summary  only 

I - data  & summary 

FLID 

11-20 

A10 

flight  id:  FLT  XYR-NN 

21-30 

A10 

date:  DD  MON  YR 

OPT 

45 

11 

0 = standard  data 

1 = Learjet  data 


averaging  interval 


(OPT  = 1 only) 


I NT 


49-50 


12 


cards  2...n  plot  request  cards 


cc 

FORMAT 

FUNCTION 

1 

11 

field  1 

3 

11 

field  2 

5-10 

16 

field  3 

15-20 

16 

field  4 

22-25 

14 

field  5 

27-30 

14 

field  6 

32-35 

14 

field  7 

37-40 

14 

field  8 

42-45 

14 

field  9 

47-50 

14 

field  10 

52-55 

14 

field  11 

57-60 

14 

field  12 

62-65 

14 

field  13 

67-70 

14 

field  14 

72-75 

14 

field  15 

77-80 

14 

field  16 

Cards  2 through  n may  appear  in  any  sequence,  however,  con- 
siderable time  is  saved  if  cards  with  same  time  limits  are 
consecutive . 

The  following  page , "Request  card  summary"  shows  the  required 
fields  for  each  plot  type. 


REQUEST  CARD  SUMMARY 
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tn 
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04 

-P 
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CL 
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X 
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AXIS  & TYPE  CODES 

The  following  axis  and  type  codes  nust  be  used  for  the  VCO 
plot  request  card. 

AXIS  CODE 


1 

Pressure 

2 

Height 

3 

Temperature 

4 

Dewpoint 

5 

LWC-JW 

6 

LWC-AUSS IE 

7 

Icing  Rate 

8 

Tacan  Bearing 

9 

Tacan  Distance 

10 

Magnetic  Heading 

11 

Velocity 

12 

Pitch 

13 

Accelerometer 

14 

LWC  Scatter 

15 

LWC  Cloud 

16 

LWC  Precip 

17 

LWC  Total 

18 

log  Z Scatter 

19 

log  Z Cloud 

20 

log  Z Precip 

21 

log  Z Total 

TYPE  CODE 

1 

2 


Scatter  Plot 
Line  Plot 
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JW-LWC  adjustment  card 

If  a VCO  plot  request  card  specifies  a LWC-JW  vs  HEIGHT 
plot,  (i.e.  h axis  code  = 5,  v axis  code  = 2,  and  type  code 
= 0)  the  next  card  must  contain  the  adjustment  parameters. 
The  required  parameters  are: 

L = number  of  levels  (ID  Maximum) 

XJ  = origin  of  the  level 
SL  = slope  of  the  level 

HT  = heig.it  (meters)  at  the  top  of  the  level 

The  card  uses  a standard  namelist  format  with  the  control 
variable  being  $ADJU3T.  If  the  option  is  not  desired  the 
card  must  be 

$ ad just  l=o , send 
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8.3  KN1UTIL  KNOLLENBERG  TAPE  UTILITY  PROGRAM 

OPERATING  INSTRUCTIONS 


CON’TF.OL  CARDS 

JQBID,CM60000,T100,TP2.  * Prob  No.  Name 

VSN , TAPEl=TAPENO , TAPE2=TAPENOIN . 

FTN,  PL:=S  99999,  SL,  R=3 , P . 

REQUEST, TAPE1 , S , HI . 

REQU EST , TA P E 2 , S , H I . 

FILE (TAP Si , RT=U , BT=K, MRL=1024 ,MBL=1024  , RB=1(  BFS=105) 
LDSET, FILES=TAPE1 , PR£SET=ZERO . 

MAP , OFF . 

LGO . 

REWIND, BPARAM, CP ARAM, DPARAM . * * 

COPY, BP ARAM . * * 

COPY, CP ARAM. ** 

COPY, DPARAM . * * 

7/8/9 

SOURCE  DECK 
7/8/9 
DATA  DECK 
7/8/9 
6/7/8/9 

*If  not  duplicating  a tape 

1.  change  TP2  to  TP1 

2.  remove  TAPE2=TAPENOIN  from  VSN  card 

3.  remove  REQUEST  TAPE 2 card 

**Renove  the  cards  if  status  word  dump  is  not  desired 
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OPTIONS 

1.  Decimal  Dump  - one  data  card 

cc  2-4  DEC 

cc  5-7  number  of  files  dumped 

cc  8-10  record  indicator 

i.e.  1 = every  record 

2 = every  other  record 

6 = every  sixth  record 

ex.  to  dump  decimally,  two  files  every  10th  record 
-DEC— 2-10  - - blank 

2.  Octal  Dump  - one  data  card 

cc  2-4  OCT 

cc  5-7  same  as  decimal  dump 

cc  8-10 

ex.  to  dump  octally,  nine  files  every  100th  record 
-OCT— 9100  - = blank 

3.  Selective  Record  Dump  - n data  cards  (Decimal  only) 
card  1 

cc  2-4  REC 

cc  5-7 (rj)  number  of  files  dumped 

cc  8-10  (rj)  record  indicator 

cards  2...n  (one  card  for  each  set  of  records) 
cc  1-6 (rj)  starting  record  to  dump 
cc  7-12 (rj)  ending  record  to  dump 

ex.  to  dump  every  10th  record 
from  125  to  350  and 
from  1000  to  1500 


-REC — 1-10 


125 350 


— 1000--1500 


- = blank 


4 . Tape  Copying  - n data  cards 


card  1 


cc  2-4 
cc  5-7 


number  of  files  to  copy 


cards  2...n  (one  card  for  each  file  to  be  copied) 

cc  1-6  (rj ) number  of  record  to  skip  before  copying 
cc  7-12  (rj)  number  of  records  to  copy 


ex.  to  create  a tape  of: 

records  701  to  750  from  file  1, 
and  records  25  to  50  from  file  3, 
the  data  deck  is  ... 


-DUP--3 
700 

24 


- = blank. 


5.  Selective  Probe  Dump  (1  data  card) 


cc  2-4 
cc  5-7 
cc  8-10 


number  of  files  dumped 

record  indicator 

i.e.  1 - every  record 

2 = every  other  record 

6 = everv  sixth  record 


cc  1 


15  selected  probe 

BL7u.'K/0  for  regular  decimal  dump  (OPTION  1) 

1 for  scatter  probe  only 

2 for  cloud  probe  only 

3 for  precip  probe  only 


A033  234 


unclassified 


DI6ITAL  PROGRAMMING  SERVICES  INC  WALTHAM  MASS  F/G  4/1 

CONTINUATION  OF  DEVELOPMENT  AND  APPLICATION  OF  OATA  PROCESSING  — ETC(U) 
JUL  76  L E BELSKY » M W FRANCIS#  F B KAPLAN  F19628-76-C-0051 

AFGL-TR-76-0182  NL 


3 3 

%33i34 


DATE 

FILMED 

2-77 


8.4  HIAClD 


KIAC1D 


The  NRI1D  program  has  been  rewritten  to  process  the  High 
Altitude  aircraft  tapes.  The  program  will  list  the  tapes  in  a 
format  similar  to  KNOLL ID . HIAClD  uses  the  standard  (210  words/ 
record)  tape  as  input.  The  program  will  produce  an  output  tape 
suitable  for  processing  with  program  RAPP  (9  words/record). 

An  editing  option  is  included  to  eliminate  certain  channels 
Actually  the  channels  specified  are  zeroed  out. 

The  median  volume  diameter  is  also  calculated  for  the  three 
probes . 


PfWIB 
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CONTROL  CARDS 


JOBNR, CM60000 , T200 , TP1 . 

REQUEST  TAPE1 , VSN= . 

REQUEST  TAPE2 , VSN= .RING. 

ATTACH  LGO, HIAC1DBINX2876 , ID=GLASS,MR=1 . 
ATTACH  A10, AlOFORMAT. 

MAP , OFF . 

LOAD, A10 . 

LGO. 

REWIND, TAPE 3. 

COPY, TAPE3 . 

EXIT, S. 

REWIND, TAPE3 . 

COPY.TAPE3. 

6/7/S/9 
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INPUT  CARDS 


Card  1 

col . 

1-6 

input 

tape 

number 

Cards  2-n 

(one 

card  per  pass) 

col . 

1-10 

start 

time 

HH: MM: SS 

11-20 

stop 

time 

HH: MM: SS 

21-25 

pass 

number 

26-  30 

averaging  interval  (second 

33 

1st  probe  to  be  edited 

34 

2nd 

•1  II 

It  II 

35 

3rd 

II  II 

•1  II 

37-38 

1st  channel 

to  be  edited 

40-41 

2nd 

II 

II  II  If 

43-  44 

3rd 

ll 

ll  II  II 

46-47 

4th 

II 

II  ll  II 

49-51 

5th 

ll 

It  II  II 

52-53 

6th 

II 

•1  II  II 

55-56 

7th 

II 

•1  ll  If 

58-59 

8 th 

II 

II  II  II 

61-62 

9th 

II 

ll  H »l 

64-65 

10th 

•1 

II  II  II 

67-68 

11th 

ll 

•1  ll  II 

70-71 

12th 

ft 

II  II  II 

73-74 

13  th 

II 

If  II  It 

76-77 

14th 

II 

II  ll  If 

79-80 


15th 


ll 


»» 
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8.5  SPANDAR  OPERATING  INSTRUCTIONS 

There  are  no  data  cards  required  to  run  this  program.  Every 
run  will  produce  an  output  tape  w ith  a tape  summary. 

CONTROL  CARDS 

JOBNR, CK65000 , T100 , TP 2 . Prob  No.  Name 

REQUEST  TAPEl , L, VSN=TAPENO . 

REQUEST  TAPE 3 , RING , VSN-LYCXXX . 

ATTACH  LGO , SPANDARBINX287  6, ID=GLASS ,MR=1 . 

FILE , TAPEl ,RT=U, BT=K,MRL=5339 ,MBL=5339 , RB=1 , BFS=536 . 


MAP , PART . 

LDSET , FILES=TAPEl , PRESET=ZERO . 
LGO. 


REWIND, TAPE8 . 
COPY,TAPE8. 


required  for  tape  listing  only 


REWIND, TAPE 2 . 
COPY , TAPE2 . 


required  for  line  printer  plot  only 


REWIND, TAPE4 . 


| required  for  punched  deck  only 


COPYCF , TAPE4 , PUNCH  ^ 


6/7/8/9 
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8 . 6 RAPP 

GLAR4,  CM75000,  T200,TP2.  2014  GLASS 

REQUEST, TAPE1,VSN=LYCXXX.  (KNOLLlD  OUTPUT  TAPE) 

REQUEST, TAPE3, VSN=LYCYYY.  (SPANDAR  OUTPUT  TAPE) 

ATTACH , LGO , RAPPBINX28  7 6 , ID=GLASS , MR=1 . 

ATTACH, CRT, OFFLINECRT, MR=1 . 

DISPOSE, TAPE39, *FL. 

LIBRARY (CRT) 

LDSET , PRESET=  ZERO . 

LGO. 

7/8/9 

DATA 

6/7/S/9 
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| 

« 

} 


CARD  1 


HEADER  CARD 


cc  1-6: 

CC  8-13: 
cc  15-24: 


cc 

25: 

cc 

31- 

40: 

cc 

41- 

50: 

cc 

55: 

cc 

56- 

65: 

CARDS 

2-N+l 

CC 

1-2: 

CC 

4-9: 

cc 

11- 

-16: 

cc 

17: 

cc 

20- 

-29: 

MISSION  ID  XYR-NN  (A6) 

A/C  TYPE  C130-E  or  C130-A  (A6) 

FLIGHT  DATE  DD0MON0YR  (A6) 

PROBE  l=SCATTER,2=CLOUD,3=PRECIP,4=TOTAL  (II) 
RADAR  OFFSET  DISTANCE  (IN  METERS) (FI 0.0) 

RADAR  CORRECTION  (IN  DB) (F10.0) 

KWAJ  (1=IF  KWAJALEIN  FLIGHT, 0=OTHERWISE) (11) 
DBl  (MINIMUM  RADAR  DETECTABLE  CONSTANT) (F10.0) 


PASS  CARDS  FOR  N PASSES 
PASS  NUBMER  (12) 

BEGINNING  PASS  TIME  IN  HHMMSS  (F6.0) 
ENDING  PASS  TIME  IN  HHMMSS  (F6.0) 
CRYSTAL  TYPE  I=ICE,  W=WATER (RAIN) (Al) 
MINIMUM  LWC  VALUE  (gm/M  )(F10.0) 


1 


1 


V 
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8.7 
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RAIN 


PROGRAM  RAIN 


CONTROL  CARDS 

JOBNM, CM70000, TP1 , T300 . 2014  GLASS 

ATTACH, LGO, RAINBINX2876 , ID=GLASS, MR=1 . 

REQUEST, TAPE1, S, VSN=TAPENO . 

ATTACH, CRT, OFFLINECRT, MR=1 . 

DISPOSE, TAPE39,*FL. 

FILE  ( TAPE1 , RT-U , BT=K , MRL=  104,  MBL=1 04,  RB=  1 , BFS“- 13) 
LIBRARY, CRT. 

LDSET , FILES=TAPE1 , PRESET=ZERO . 

LGO. 

7/8/9 

DATA  CARDS 
6/7/8/9 


- 


DATA  CARDS 


CARD  1 

cc  1-2: 
cc  4-13: 


CARD  2- (F+l) 

cc  1-6: 
cc  10-15: 
cc  20-25: 


HEADER  CARD 

GROUPING  INTERVAL  (SECONDS)  (USUALLY  10) 
DATE  DD0MON0YR  @ = blank 


F = NUMBER  OF  FILES  ON  TAPE 

FILE  START  TIME  (OCTAL)  HHMMSS 
(K)  NUMBER  OF  SIGNIFICANT  BITS  (15) (r j ) 
'F)  RANGE  SWITCH  SETTING  OF  PEAK  LOAD 
VOLTMETER  (rj) 


